Master crossover functions for the one-component fluid "subclass" 
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Introducing three well-defined dimensionless numbers, we establish the link between the scale 
dilatation method able to estimate master (i.e. unique) singular behaviors of the one-component 
fluid "subclass" and the universal crossover functions recently estimated [Garrabos and Bervillier, 
Phys. Rev. E 74, 021113 (2006)] from the bounded results of the massive renormalization scheme 
applied to the (n)-model of scalar order parameter (n = 1) and three dimensions (d — 3), 
representative of the Ising-like universality class. The master (i.e. rescaled) crossover functions are 
then able to fit the singular behaviors of any one-component fluid without adjustable parameter, 
only using one critical energy scale factor, one critical length scale factor, and two dimensionless 
asymptotical scale factors, which characterize the fluid critical interaction cell at its liquid-gas critical 
point. An additional adjustable parameter accounts for quantum effects in light fluids at the critical 
temperature. The effective extension of the thermal field range along the critical isochore where the 
master crossover functions seems to be valid corresponds to a correlation length greater than three 
times the effective range of the microscopic short-range molecular interaction. 

PACS numbers: 64.60.Ak., 05.10.Cc, 05.70.Jk, 65.20. +w 



1. INTRODUCTION 

The universal features of three-dimensional (3D) Ising- 
like systems are now well-established by the renormal- 
ization group approach (author?) [l|] of the classical-to- 
critical crossover behavior (author?) Q. In this theo- 
retical context, it was possible to estimate the complete 
functions which interpolate between the critical behavior 
(controlled by the non-trivial (Wilson-Fisher) fixed point 
(author?) 0,11) and a classical behavior (controlled by 
the Gaussian fixed point) . Such interpolating theoretical 
expressions were customarily named classical-to-critical 
crossover functions. The corresponding crossover within 
the critical domain was referred as the critical crossover 
author?) 0,0], or the asymptotic crossover (author?) 



Our present interest is restricted to the dimension- 
less expressions derived from a massive renormalization 
(MR) scheme applied to the §\ (n) model, for three- 
dimensional systems (d = 3) and scalar order parame- 
ter (n = 1) (author?) @, @, 0, E3] (d and n are the 
dimensions of the space and order parameter density, 
respectively, which characterize each universality class 
(author?) Q). In that specific renormalization group 
approach (author?) the Ising-like universality is 

linked to the existence of a unique non-trivial fixed point. 
For convenient simplification in the following presenta- 



tion, the complete Ising-like universality class is labeled 
{$3 (l)}-class (with reference to the $^ =3 (n = l)-model), 
while the one-component fluid "subclass" made of all one- 
component fluids is labeled {l/}-subclass, with the ob- 
vious relation {l/}-subclass C {$3 (l)j--class. 

The introduction of the system-dependent parame- 
ters for the practical use of the theoretical functions 
was discussed in a detailed manner in Refs. (author?) 
0, [l(| E2, G3 • More generally, the dimensionless forms of 
the theoretical expressions must be used to fit the exper- 
imental results in order to preserve both the number and 
the critical scaling nature of the fluid-dependent factors 
which are free in the massive renormalization scheme. 
Indeed, it was precisely shown in Ref. (author?) 
hereafter labeled I, that the Ising-like universal features 
(author?) [II, estimated in the Ising-like preasymptotic 
domain close to the non-trivial fixed point, require to 
characterize each system along its critical isochore (the 
thermodynamic line equivalent to ft. = 0) using four pa- 
rameters. 

Two of them are dimensional parameters, namely the 
critical temperature T c which acts as energy unit (in- 
troducing the universal Boltzmann constant ks) to ex- 
press the Hamiltonian in dimensionless form, and the 
unknown inverse coupling constant (go) of the fourth- 
order term of the dimensionless Hamiltonian. As a mat- 
ter of fact, any Hamiltonian representative of a physical 
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system at near criticality, such as a one-component fluid 
near its liquid- vapor critical point, is driven to the non- 
trivial fixed point under the action of the renormaliza- 
tion transformations (author?) Q. Due to the fact that 
renormalizable field theories are short-distance insensi- 
tive, universality emerges in a regime £A 3> 1 in which 
the correlation length £ is much larger than the micro- 
scopic scale, which plays the role of the inverse wavenum- 
ber cut-off (Ao) in the renormalization scheme. This 
universality is non-mean field like in nature (at least for 
the three-dimensional systems which are of present in- 
terest), because the actual molecular interaction range 
at the microscopic scale of the physical system cannot be 
completely eliminated. (Ao) - remains the single natu- 
ral length unit in the theoretical scheme. Indeed, (go) -1 
takes the convenient length dimension at d — 3 to act 
as adjustable length unit, and to express the correlation 
length in dimensionless form. 

The other two parameters are dimensionless coeffi- 
cients, namely the scale factors d and tp, which provide 
the analytical (linear) proportionality between the two 
dimensionless physical fields At* and Ah* of the Ising- 
like fluid and the two renormalized relevant fields t and h 
of the <5^_ 3 (n = l)-model, respectively (author?) 0,0] 
(using customarily field notations, see below). In such a 
situation, the universal features close to the non-trivial 
fixed point are estimated in conformity with the so-called 
two-scale-factor universality, where only two asymptotic 
critical exponents and one confluent exponent are inde- 
pendent. Then, the lowest value A ~ 0.51 (author?) 
[14 ] of the confluent exponent characterizes the correc- 
tions to scaling due to one possible irrelevant field (au- 
thor?) [15|). In order to maintain the coherence with 
the previous presentation of these universal features given 
in Refs. (author?) [T3, 0], we have here also selected 
v ~ 0.630 and 7 ~ 1.240 (author?) P3] as indepen- 
dent leading exponents attached to the correlation length 
^th (t)) and the susceptibily Xth (*) along the critical iso- 
chore (h = 0), respectively. 

The finite scale (Ao) 1 is generally unknown for a real 
microscopic interaction at short range distance in pure 
fluids. Simultaneously the macroscopic size L of the 
fluid sample should be larger than £, i.e. L £ — > 00, 
then a special attention to account for extensive nature 
of the thermodynamic properties of the physical system 
is needed. Moreover, thanks to the general point of view 
of the thermodynamics for 3-D systems, the dimension- 
less forms of any physical density variable ^ (where X 
is the total extensive variable and V oc L d is the total 
volume of the system) can be obtained without reference 
to the unknown wavelength number go defined at the 
critical point. Indeed, introducing the total amount of 
matter N [or the total mass M = Nm p , where nip is the 
mass of the particle, while the subcript p refers to a par- 
ticle property] of a system filling a total volume V, the 
dimensionless order parameter conjugated to the dimen- 
sionless ordering field can always be defined if the amount 
of matter in a reference volume is known (here the refer- 



ence volume can be chosen for example as the volume of 
a mole, a particle, a cell lattice, a mass unit of matter, 
etc.). Therefore, any reference length a , defined such as 
no is the amount of matter in the volume {ao) d , can be 
used as explicit length unit for the thermodynamic and 
correlation functions. Thus the massive renormalization 
scheme generates a third adjustable dimensionless scale 
factor - namely Uq — <?o a o - which relates the dimen- 
sionless correlation length ^- of the physical system to 
the corresponding theoretical function derived from the 
massive renormalization scheme. As a correlative result, 
when ao takes its physical sense to represent the effective 
range of the microscopic molecular interaction between 
the no particles, i.e. ao oc (Ao) while no oc coordina- 
tion number, the Ising-like singular nature of the physical 
system can be characterized by a set of three dimension- 
less scale factors ip}. However, in such a three- 
parameter characterization of the physical system, it is 
then essential to recall that the theoretical estimations 
of the universal features are only valid within the Ising- 
like preasymptotic domain. In this preasymtotic domain, 
each dimensionless theoretical function can be approxi- 
mated by its restricted asymptotic form as a two-term 
Wegner-like expansion, leading to three independent crit- 
ical exponents (i. e. our selected set {v, 7, A} in present 
work) . 

Indeed, in the seventies, it was clearly shown by ex- 
perimentalists that the singular properties of pure flu- 
ids close to their liquid-gas critical point were satisfied 
by power laws with universal features comparable to the 
ones estimated for the uniaxial three-dimensional Ising 
system used as a predictive model (for a review see for 
example (author?) [lfij|). It was then revealed that two 
independent leading amplitudes, attached to the univer- 
sal values of two independent critical exponents, are the 
only two fluid-dependent parameters necessary for char- 
acterizing the asymptotic singular behavior of each one- 
component fluid. Therefore, selecting the dimension- 
less correlation length £* (At*) and the dimensionless 
isothermal compressibility K* T (At*) in the homogeneous 
domain (At* > 0) along the critical isochore (Ap* = 0), 
each Ising-like critical fluid can then be characterized by 
the related leading amplitudes £ + and T + (using stan- 
dard notations for critical fluids (author?) [13]). This 
asymptotic situation characterized by two dimensionless 
leading amplitudes £ + and T + was in conformity with 
the two-scale-factor universality expected for all systems, 
with short-ranged interaction, and which have an isolated 
transition point. 

Correlatively, it was demonstrated (author?) [l^] that 
the two-scale-factor universality related to the Ising- 
like nature of the critical phenomena in pure fluids are 
"observed" in a very limited range of temperature and 
densities around their liquid-gas critical point. Obvi- 
ously, since the asymptotical critical domain associated 
to this limit is so narrow that experiments are diffi- 
cult to achieve, it was fundamental to account for the 
possible nonuniversal character of the system through 
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the confluent singularities in the corrections to scaling 
(author?) (ignoring here the background contri- 
butions which are significative only in the case of spe- 
cific heat (author?) [12]). That leads to express the 
singular properties as truncated forms of the Wegner- 
like series. That precisely corresponds to the Ising-like 
limit of the asymptotic crossover mentionned just above 
and investigated in details in the renormalization theory 
where the resummation of the Wegner-like expansions 
should yield complete crossover functions from asymp- 
totic (Ising-like) critical behavior to the classical (mean- 
field) critical behavior. Different theoretical approaches 
have been adopted by many investigators to obtain ex- 
plicit solutions resumming the complete Wegner series 
(see for exemple a review in Ref. (author?) fff] for their 
application to the fluid case). The practical results essen- 
tially depend on the approximations used in the renor- 
malization scheme and the way to account for the cutoff 
effects. Despite these technical differences to treat the 
asymptotic crossover, the Ising-like universal feature was 
related to the lowest confluent exponent A where only 
one (fluid-dependent) confluent amplitude is needed to 
characterize the first order term of the confluent correc- 
tion to scaling (author?) [l5j]. In the characterization of 
each critical fluid where the leading amplitudes £ + and 
T+ are selected in conformity with the asymptotic two- 
scale-factor universality, then it was necessary to add the 
confluent amplitude a+ of the first order correction term 
for the susceptibility (related to the confluent amplitude 
a~£ of the correlation length by the universal value of the 

ratio -qr « 0.68 (author?) [14]). The resulting ampli- 

a x 

tude set {£ + ,r+,a+} defines the complete asymptotic 
crossover of each one-component fluid. This amplitude 
set is Ising-like equivalent (in quantity and nature) to the 
previous scale factor set ip} used to characterize 

the Ising-like critical behavior whitin the preasymptotic 
domain of each one-component fluid. 

As recalled in I, this three-parameter description of the 
asymptotic singular behavior of the correlation length 
and the susceptibility of xenon was studied (author?) 

using the crossover functions initially derived by Bag- 
mils and Bervillier (author?) @, Q from the massive 
renormalization scheme. In this pioneering study of the 
crossover, for the first time the minimal quantity of Ising- 
like non-universal parameters of xenon was introduced as 
a set of a single wavelength (defined at the critical point), 
and two dimensionless scale factors expressing the ana- 
lytical approximation between the two relevant scaling 
(thermal-like t and magnetic-like h) fields and the corre- 
sponding physical (At* and Ah*) fields (author?) 
Subsequently, theoretical and numerical approaches ap- 
plied to the asymptotic crossover description of the singu- 
lar behavior observed in pure fluids, have confirmed this 
characterization with thr ee p arameters (see for example 
Refs. (author?) 0, E3, M, EI 0, El and a review in 
Ref. (author?) |<f)7 

Today, with the appropriate introduction of two fluid- 



dependent factors in conformity with the two-scale-factor 
universality, any theoretical function which fits the tem- 
perature dependence of the effective critical exponent 
along the critical isochore may be made universal by 
simply rescaling the temperature distance to the critical 
temperature, as initially proposed by Kouvel and Fisher 
(author?) j25| who introduced a single crossover tem- 
perature scale At£-. Unfortunately, the unsolved prob- 
lems in these theoretical approaches remain the validity 
of the linear approximations of two relevant fields (which 
correctly introduce the two system-dependent scale fac- 
tors) , the importance of the neglected analytical and non- 
analytical corrections, and, more generally, the estima- 
tion of the extension range in temperature and densities 
around the liquid gas critical point where the Ising-like 
universal features should be observed. 

To complete the above introduction of the non- 
universal character of the asymptotic crossover in pure 
fluids, we also recall that, at the beginning of the eighties, 
the description of the behavior of the singular thermody- 
namic properties at finite distance from the liquid gas 
critical point was also made using the theoretical formu- 
lation of the nonasymptotic crossover from a regime of 
Ising-like scaled behavior to another regime in which the 
critical anomalies due to large fluctuations are ignored 
(author?) 0, HI, EE HI, The common attempt 

to address this problem was based on the classical-to- 
critical crossover description of the free energy density. 
Indeed, this approach is useful for better understand- 
ing of crossover critical phenomena in "complex" fluids 
where the character of the crossover reflects an interplay 
between Ising-like universality caused by long-range fluc- 
tuations and a specific supramolecular structure charac- 
terized by an additional nanoscopic or mesoscopic length 
scale (which can then differ significantly from (Ao) ). 
Therefore, while the Ising-like two-scale-factor universal- 
ity was similarly accounted for introducing the two di- 
mensionless parameters of proportionality between the 
respective relevant (physical and renormalized) fields (for 
example Ct (~ i?) and c p (~ V) m the notations of Refs. 
(author?) [H, H3]), the fundamental difference with 
the asymptotic crossover description comes from the in- 
troduction of two independent dimensionless parameters 
(for example u and A in the notations of Ref. (author?) 
[23]), m order to control this nonasymptotic crossover 
character in complex fluids. However, in an application 
related to the pure fluid case, it is not necessary to in- 
troduce an additional mesoscopic length scale to account 
for the realistic microscopic situation in one-component 
fluids (author?) [3l|. Such pure fluids can then be 
assimilated to Lennard- Jones-like fluids when they are 
made of atoms or highly centro-symetrical molecules, or 
to short-distance associating fluids when they include 
more sophisticated short-range molecular interactions be- 
tween unsymetrical molecules, polar molecules, bonding- 
like molecules, etc. Moreover, the representation of the 
experimental phase surface of any pure fluid by a van der 
Waals-like equation of state is not accurate either close 
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or far away from the critical point (since the van der 
Waals equation of state is theoretically justified only for 
infinite range of the molecular interaction) . As a final re- 
sult, the fluid-dependent parameters needed to describe 
the classical behavior of the free energy density have no 
quantitative signification. The nonuniversal complexity 
of the pure fluid was then accounted for by introducing a 
significant number of adjustable parameters whose cou- 
pling with the two dimensionless crossover parameters 
u and A can not be completely defined. Therefore, in 
spite of the correct introduction of a crossover function in 
the definition of variables and thermodynamic potentials, 
the only founded theoretical challenge of the nonasymp- 
totic crossover applied to the one-component fluids re- 
mains to account for the correct Ising-like universal fea- 
tures with a single crossover scale. The uniqueness of the 
crossover scale can thus be defined introducin g an arbi- 
trary fixed value of the product uA (author?) [30| . The 
set {u (or A) , c t ,c p } appears "Ising-like" equivalent to the 
set {uq, i9, ip}- This nonasymptotic crossover (which thus 
must match the asymptotic critical crossover close to the 
Wilson-Fisher fixed point (author?) [32]) has to be not 
completely solved in regards to the most recent theoreti- 
cal predictions of universal exponents (author?) fl4H33T| 
and universal amplitude ratios (author?) [jiKlJ]. More- 
over, the introduction of the single crossover parameter, 
which is then related to the mean-field concept of the 
Ginzburg number (author?) (28|, add conceptual diffi- 
culties to understand the role of real microscopic param- 
eters controlling a true rescaled universal behavior in the 
whole crossover region (author?) 0,0, HI- 

Finally, since the van der Waals dissertation, the real 
difficulty for scientists interested in liquid-gas critical 
phenomena in pure fluids, comes from the nonclassical 
(i.e. renormalizable) theories which are not able to pre- 
dict the location of the critical point, while the classical 
theories provide its uncorrect location. Such a difficulty 
has generated a crucial experimental challenge where the 
determination of the two characteristic leading ampli- 
tudes and the characteristic crossover parameter of each 
pure fluid, and alternatively but equivalently, the local- 
ization of its liquid-gas critical point on the p, V, T phase 
surface, remain mandatory. 

Based on this recurrent situation, an alternative phe- 
nomenological way to characterize the asymptotic singu- 
lar behavior of the one-component fluids was also for- 
mulated by Garrabos (author?) [33] as follows: "If you 
are able to locate a single liquid gas critical point on 
the experimental p, Vp, T phase surface of a fluid parti- 
cle of mass m p , then you are also able to describe the 
asymptotic crossover around this isolated point", p is 
the pressure, T is the temperature, Vp = ^ = — is 
the volume of the particle, and p is the (mass) density. 
Accordingly, a minimal set Q™™. made of four critical 
coordinates (author?) [35] [see below Eq. ([!])], pro- 
vides unequivocal determination of four (two dimensional 
and two dimensionless) scale factors [see below Eqs. ([3]) 
to l[7])]. Then a scale dilatation method of the physi- 



cal fields can be used to observe and quantify the mas- 
ter (i.e., unique) asymptotic crossover behavior of the 
{f/}-subclass (author?) [H, The two dimensional 
critical parameters, noted (f3 c ) 1 and a c , take appropri- 
ate energy and length dimensions, respectively to reduce 
the physical variables, the thermodynamic functions, and 
the correlation functions. The two dimensionless critical 
numbers, noted Y c and Z c , are well-defined characteris- 
tic parameters of the critical interaction cell of volume 
{ct c ) d ■ An additional adjustable parameter, noted A* e , 
accounts for quantum effects in light fluids at the criti- 
cal temperature (author?) [38]. Conversely, when Q™™ 
and A* e were known for the selected fluid, the asymptotic 
master behavior characterized by three master (i.e. con- 
stant) amplitudes was used to calculate the amplitude set 
{£ + , r + , a+ } which characterizes the asymptotic singular 
behavior of this fluid. In addition to this intrinsic pre- 
dictive power, another important characteristic attached 
to the scale dilatation method was the Ising-like anal- 
ogy in its formal introduction of the two dimensionless 
scale factors Y c and Z c and the corresponding ones i? 
and ip introduced by linear approximations in the mas- 
sive renormalization scheme. 

As a matter of fact, for each selected fluid belonging 
to the {l/}-subclass, this analogy can be useful to pro- 
vide explicit estimation of the unknown scale factor set 

{(ffo) 1 ; V'l [or {u^, ?/>}] of the theoretical crossover- 
functions (using then, the thermodynamic length scale 
unit a c of the selected one-component fluid as a reference 
length ao). Especially in the case of the unique form of 
the mean theoretical functions estimated in I (which in- 
corporates the error-bar propagation of the min and max 
crossover functions revisited in (author?) [13]), we can 
formulate the unambiguous modifications of the theoret- 
ical crossover functions for the {$3 (l)j--class to exactly 
match the master two-term Wegner-like expansions valid 
within the Ising-like preasymptotic domain of the {1/} 
-subclass. 

These formulations were used to study the correla- 
tion length in the homogeneous domain of seven one- 
component fluids (author?) [3^| and the squared cap- 
illary length in the non-homogeneous domain of twenty 
one-component fluids (author?) [13] ■ Similarly, a re- 
cent application to the practical parachor correlations 
(i.e., equations expressing surface tension as a power law 
of the density difference between coexisting gas and liq- 
uid phases), have shown that the corresponding master 
form acts as a universal equation of state for the inter- 
facial properties (author?) [4lJ. Now, our present ob- 
jective is to achieve the complete uniquevocal link be- 
tween these updated results of I and the scale dilatation 
method to predict the master singular behavior of the 
{l/}-subclass. For these studies, the analytical relations 
between the relevant scaling fields of both descriptions 
must be defined. 

The paper is organized as follows. In Section 2 the 
master description of the universal features within the 
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Ising-like preasymptotic domain is recalled. First, start- 
ing from the four critical coordinates of the critical point, 
we define four scale factors which are needed to un- 
ambiguously determine three dimensionless amplitudes 
which characterize the Ising-like preasymptotic domain of 
each one-component fluid. Second, we show the master 
singular behavior of the isothermal compressibility, ap- 
plying the scale dilatation method to the related physical 
quantities. That complete the master sigular behavior of 
the correlation length in conformity with the two-scale- 
factor universality of the {$3 (l)}-universality class. In 
Section 3, a brief presentation of the theoretical crossover 
functions for the correlation length and the susceptibil- 
ity in the homogeneous phase is given to demonstrate 
the analytical matching with the master singular behav- 
ior provided by the scale dilatation method. Introducing 
three well-defined dimensionless numbers characterizing 
the {l/}-subclass, the unequivocal link between three 
theoretical amplitudes, which characterize the {$3 (1)}- 
universality class, and three master amplitudes, which 
characterize the {l/}-subclass, is given before concluding 
in Section 4. Two appendices deal with first, the equiva- 
lence between different one-parameter crossover models, 
and second, the determination of the crossover parameter 
beyond the preasymptotic domain using the well-known 
linear model of the parametric equation of state with ef- 
fective exponents. 



oir.c = |^7 C - The subscript c refers to a critical quantity. 



From Eq. ([!]), we can construct a more convenient set, 

Qt*={{0 c T\a c ,Y c ,Z c ) (3) 
making use of the following four scale factors 



(/3 C ) 1 = k B T c ~ [energy] 



O'e 



knT c 



p c n 



kBT c 



[length] , 



[dimensionless] 
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[dimensionless] 



(4) 



(5) 



(6) 



(7) 



(Pe) and a c are used to express dimensionless quanti- 
ties. a c is a measure of the effective range of the micro- 
scopic short-range molecular interaction (Lennard- Jones 
like in nature) (author?) [3l[. Z c is the critical com- 
pression factor, while Y c — or, c — 1. In the above dimen- 
sionless form of the thermodynamic functions normalized 
per particle, 4- is the number of particles in the volume 



2. MASTER SINGULAR DESCRIPTION OF 
THE ONE-COMPONENT FLUID SUBCLASS 

2.1. The minimal set of critical parameters 

For the {l/}-subclass, it was hypothesized (author?) 
[35I ] (author?) [36j that all the information needed to 
characterize non-quantum fluid critical phenomena is 
contained within the four critical parameters needed to 
localize the single critical point and its tangent plane on 
the experimental phase surface of normalized equation of 
state (p, Vp, T) = (the needed supplementary in- 
formation to characterize quantum fluids is given in Ref. 
(author?) Q; see also below Eqs. © and (|T0]l), This 
minimal set of four coordinates reads as follows 



5 c.a« 



\T c ,p c ,vp, c ,~f' e \ 



(1) 



where 



Vs,c = -*r = is the critical volume per particle 
(V is the total volume, N c is the total critical number of 
particles, and p c is the critical density), and 



7c = ^ 



dp 
df 



C ;T=T C 



dp s 



(IT 



(2) 



T=T C 



is the common critical direction of the critical isochore 
and the saturation pressure curve at the critical point, in 
the p; T diagram. 7 C is related to the Riedel factor (au- 



v c ,i = {a c ) c 



(8) 



which corresponds to the volume of the critical interac- 
tion cell (author?) [35| . 

This actual set Q™ n (made from measured critical pa- 
rameters), refers to the characteristic range of the mi- 
croscopic molecular interaction in "classical" (i.e. non- 
quantum) fluids [here the molecular interaction range is 
measured by a c of Eq. {HJ)]. To include quantum fluids 
in the one-component fluid subclass (author?) [38|, we 
need the phenomenological introduction of a supplemen- 
tary adjustable parameter, noted A* e , which accounts for 
the quantum effects at this microscopic length scale of 
the effective molecular interaction. The (dimensionless) 
parameter A* e (author?) [38] is given by 



A* 



1 + A C 



with 



A, = A 



A 



T.c 



a. 



(9) 



(10) 



thor?) [4J], a R ,, 



dlogT 



X q j (with Xqj > 0), is thus a non universal adjustable 
number which accounts for statistical contribution due to 
the nature (boson, fermion, etc.) of the quantum particle. 
At c = — r is the de Broglie thermal wave- vector 

at T — T c , hp is the Planck constant (the subscript P 
is here added to make a distinction with the theoretical 



T=T C 



through the relation ordering field noted h) . 
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2.2. 



Thermodynamic characterization of the 
critical interaction cell 



Nrn^ 



We introduce the (mass) density variable p — — 
j£ and we consider the usual compression factor 



Z 



pV 



pm ¥ 



Nk B T pk B T 



(11) 



generally expressed in thermodynamic textbooks (au- 
thor?) [Ill as a function Z(T*,p) of the two dimen- 
sionless variables T* = ip- and p = — . Here we note 

Ic r Pa 

the distinction underlined using superscript asterisk for 
a dimensionless quantity obtained only from(/3 c )~ 1 and 
a c units, and decorated tilde for a dimensionless quan- 
tity which can refer to a specific amount of matter, then 
introducing also the critical density p c . Practically, the 
two dimensionless critical parameters 



jj.c 



dZ_ 

df* 



P=Pc 



J CP 



dZ' 

dp / T » =T , 



= -Z r 



(12) 



(13) 



J CP 



are the two preferred directions (author?) [44| of the 
characteristic surface related to the total Grand poten- 
tial J (T,V, pp), expressed per particle. p p is the chemi- 
cal potential per particle related to the specific (i.e., per 
mass unit) chemical potential p p by p p = ^ (where the 

subscript p refers to a specific property). Therefore, it is 
essential to note that y* p c — Y C Z C and Zp c = —Z c are the 
dimensionless forms of two characteristic molecular (i.e., 
per particle) quantities. 

As a matter of fact, when we consider the ther- 
modynamic description of a one-component fluid at 
constant volume of matter, the total Grand poten- 
tial J (T, V, pp) = —p (T, pp) V takes, alternatively but 
equivalently, the role of the total Gibbs free energy 
G(T,p,N) = pp(T,p)N usually considered in the 
thermodynamic description of a one-component fluid 
of constant amount of matter. The external pressure 
p (T, pp) = of the container maintained at constant 
volume, in contact with a particle reservoir, is then 
the thermodynamic potential equivalent to the molec- 
ular chemical potential p p (T, p) = % of the fluid main- 
tained at constant amount of matter, in contact with a 
volume reservoir. Therefore, considering the normaliza- 
tion per particle of the thermodynamic description of a 
one component fluid at constant volume, the molecular 
(i.e., per particle) Grand potential reads, jp tVp =cte. (T) — 
—p(T,pp)vp. Using the associated opposite Massieu 

form, Zp^ Vp=cte . = — ^ JpiT, ^r cte ' , and the "universal" 

Boltzmann constant ks as unique unit, we obtain the 
following dimensionless form 



D,tip=cte. 



p (r, pp) v p 

k B T 



(14) 



which demonstrates that the compression factor Z of a 
constant amount of fluid matter maintained at constant 
volume (i.e. a one-component fluid monitored by the 
temperature along an isochore) is indeed a dimension- 
less molecular potential (author?) [45]. For the criti- 
cal filling N — N c of this isochoric container, we obtain 

* * _ r (rTnl i 

Z c . 



( Jp,Up=cte. \ 



N=N C 



V 



Pr 



P{T) 
T 



Here, 



P(T) 
T 



P=Pc 



P=Pc 



acts as first characteristic (i.e., inde- 



pendent) equation of state for a critical isochoric fluid, 
where the two extensive variables V and N c are fixed 
[i.e., a critical fluid at p = 1 in contact with a thermo- 
stat (i.e. an energy reservoir) of constant energy ksT]. 
Multiplying the particle property y% c by the number of 
particle in the critical interaction cell, it appears that 



the critical quantity Y c 



dT* 



is read- 



cp 



ily a characteristic parameter of the critical interaction 
cell. 

Now considering a critical isothermal fluid where the 
two variables V and T c are fixed (i.e., a critical fluid at 
T* = 1, filling a constant total volume thermostated at 
constant critical energy fcsT c , in contact with a particle- 
obtain - ( k^ssil )* = z 7 



reservoir 



we 



l "\ 



Ip(pp)] t = 



T=T a 



L v e . Here 



p{pp) 



acts as sec- 

T=T C 

independent) equation of state 



ond characteristic (i.e. 
for a critical isothermal one component fluid. In such a 
thermostated container at fixed total volume, we under- 
line the fact that the only independent extensive variable 
to monitor the thermodynamic fluid state is the number 
of particles TV which fixes the equilibrium mean value of 
the molecular chemical potential pp. For N = N c , at 
T* = 1 (i.e. the critical point condition), the critical 
chemical potential per particle takes the value pp. c , such 



that 



T=T C 



Z c . Within the 



critical interaction cell filled with -g- particles, the nor- 
malized Grand potential takes the master critical value 

± U \ =1 

Therefore, as an essential microscopic meaning related 
to Eq. |(8]), we note that the critical set Q™ m of Eq. 

, characterizes the master thermodynamic information 
contained in the critical interaction cell volume of each 
one-component fluid at the critical point. 

Finally, we summarize the two main constraints for 
the thermodynamic description of a one-component fluid 
near its gas- liquid critical point: 

i) The dimensionless reduction of the variables is 
mandatorily made by using the two dimensional factors 
{(i c ) 1 and a c of Eq s, Q and ([5]), respectively (see also 
Ref. (author?) pj); 

ii) The thermodynamic properties expressed per par- 
ticle are better suited to understand the microscopic na- 
ture of the two dimensionless numbers Y c and Z c . That 
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leads to express dimensionless properties from reference 
to the ones estimated for the volume of the critical inter- 
action cell. Then the thermodynamic origin of the dimen- 
sionless master (i.e., unique) constants is well-identified. 



2.3. The relevant physical fields crossing the 
liquid-gas critical point 

Such a constrained dimensionless thermodynamic de- 
scription is appropriately obtained from the Grand 
canonical statistical distribution, considering a one- 
component fluid in contact with a "particle-energy" reser- 
voir maintained at constant total volume V, Selecting 
the thermodynamic nature (fixing, either the energy level 
/csT, or the particle amount N) of the reservoir to reach 
the critical point (either at constant critical density, or 
constant critical temperature) , the normalized thermody- 
namic potential is then related to the intensive quantities 

In addition to the tempera- 



p(T) 

T 



or 



P=p c 



T 



T=T C 



ture variable conjugated to the total entropy, the other 
natural (intensive) variable is the chemical potential per 
particle fip, conjugated to the natural fluctuating total 
number of particles N (leading to the fluctuating num- 
ber density n = ^). Therefore, the two relevant physical 
fields, either to express the finite distance to the critical 
point, or to cross it, along the critical isochore and along 
the critical isotherm, are 



At* = k B f3 c (T - T c ) 



(15) 



and 



Ah* = /3 C (fip — fj,p iC ) (16) 

respectively. Using the thermodynamic description per 
particle, the order parameter density is then proportional 
to the critical number density difference n — n c (n c = ^f- 
is the number density), and the associated dimensionless 
order parameter density is given by (author?) [aEE]: 



Am* = (n - n c ) (a c ) c 



(17) 



We retrieve the distinction (using superscript asterisk or 
decorated tilde), either between Ah* [see Eq. lfl6j) ]. and 



A/2 = {p p - /i p , c ) 



(18) 



which show that the dimensionless isothermal suscepti- 



bilities Xt 



a(Am') 



a factor 



d(Ah 

Equations 



and xt 



d(Ap) 
9(Ap) 



differ by 

f 



2Q| and ([2T|) illustrate the 



primary role of Z c in the dimensionless form of thermo- 
dynamics, due to the fact that (Z c ) , i.e., the parti- 
cle number within the critical interaction cell volume, ac- 
counts for extensivity of the critical fluid. 



2.4. The scale dilatation method for the 
{l/}-subclass 

A detailed presentation of the scale dilatation method 
can be found in references (author?) [H, HI, EE H 5*1 . 
Hereafter we only recall the main features which close the 
master description of the singular behaviors of the {1/}- 
subclass within the preasymptotic domain (with 7, v, 
and A selected as independent critical exponents). The 
scale dilatation method uses explicit analytical transfor- 
mations of each physical field At* and Ah* given by the 
equations 



T q * f = T* = Y c I At* 



-T_/* 



(A* qe ) 2 H* = (A* qe ) 2 (Z c )-?\Ah* 



(22) 



(23) 



where T* ( = T* is the renormalized thermal field, and H* [{ 
is the renormalized ordering field. The subscript qf dis- 
tinguishes between a quantity which refers to a quantum 
fluid (i.e., A* e ^ 1) from the one which refers to a non- 
quantum fluid (i.e., A* e = 1) (author?) [3S]. Accord- 
ingly, the analytic transformation between the physical 
order parameter density Am* and the renormalized or- 
der parameter density reads as follows (author?) 

SHE! 



M* q( = A* qe M* =A* qe (Z c y 



I Am* 



(24) 



Introducing then the dimensionless correlation length 
C* = q , the renormalized correlation length £* f is given 
by the equation 



or between Am* [see Eq. Q17p]. and 



Ap = 



P- Pc 



(19) 



where A// and Ap were customarily defined in a critical 
fluid description using specific properties and practical 
dimensionless variables x = (see, for example, Refs. 
(author?) @, [lB|). The corresponding relations can be 
expressed as follows, 



Ah* 
Am* 



Z c Ap 
^P 



(20) 
(21) 



C=(A. 



qe) 



(25) 



which preserves the same length unit for thermodynamic 
and correlations functions (with i* = £* for the non- 
quantum fluid case). 

The master asymptotic singular behavior of £* f (T*) 
was studied in (author?) [39]. Specifically, the observed 
asymptotic divergence of £* [t was represented by the fol- 
lowing (two-term) Wegner expansion 



f* _ 7+ 



(r * r 



rl,+ 



A 



(26) 
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Z ± 


{^ + } 


Z x± 




pO,± 


P i,± 




Z+ = 0.570481 


ZJ = 0.291062 




Z X A = 0.37695 
Z^~ = 0.37695 




? ±=« c? ± = a c A* e (y c )-" 2± 


a±=4'±(y c ) A 




Z+ = 0.119* 


Z~ = 0.0248465 


= 0.555* 


Z)l~ = 2.58741 




r± = (A* e ) d - 2 (z c )-i(^c)-^± 


a±=^±(y c ) A 






Z% = 0.105658 
Z~ = 0.196829 




Z]j + = 0.522743 
z)f = 0.384936 






a±=^±(y c ) A 






Z M = 0.468* 




Z X M = 0.4995 


&P*LV 


B=(k* qe )- 1 {Z c yh(Y c )PZ M 


a M = zi x (y c ) A 



Table I: Master (Vqi) and physical (P ) amplitude values of the singular behavior of the correlation length (lines 2 and 3), the 
susceptibility (lines 4 and 5), the specific heat (lines 6 and 7) and the order parameter density (line 8), along the critical isochore 
of any one-component fluid; columns 2 and 4: independent master amplitudes [see Eq. (|40 j) : columns 3 and 5: amplitude values 
in conformity with the theoretical universal features estimated within the Ising-like preasymptotic domain (author?) [Tol. [T^|: 
asterisk indicate the "experimental" master values estimated using xenon as a standard critical fluid (see Refs. (author?) 
[36L l37l. l48l|);columns 7 and 8; corresponding physical amplitudes when Q™ n = {(/3 C ) _1 , a c ,Y c , Z c } and A*, e are known for the 
selected one-component fluid. 



where v = 0.6303875 and A = 0.50189 (author?) [II]. 
The leading amplitude Z^ = 0.570481 and the first 

confluent amplitude Z^' + = 0.37695 have master (i.e. 
unique) values for the {l/}-subclass. The associated 
asymptotic singular behavior of the physical correlation 
length was given by 



£ex P (At*) = £ + (At*)-" 1 + 4 (At 



(27) 



Therefore, the term to term comparison of (master) Eq. 
(|26|) and (physical) Eq. |27j) . results in the following 
amplitude combinations 



Ot c 



- Z 1 ' 



(28) 



(29) 



Applying now the scale dilatation method to any phys- 
ical (thermodynamic) property P(Ar*), the master sin- 
gular behavior for the renormalized (thermodynamic) 
property V* { (T*) can be also observed and represented 
within the preasymptotic domain by the restricted ex- 
pansion 



4'± (T*) A 



(30) 



where Zp and Zp^ are two master constants for any 
one-component fluid (see Table [J. To close the mas- 
ter description in conformity with the universal features 
estimated within this Ising-like preasymptotic domain, 
we complete the representation of the master correlation 
length with the one of the master susceptibility X* { ob- 
tained from master order parameter density Al* f , and 
master ordering field 7i* f , using the thermodynamic def- 
inition, X* { 



X* { is related to the dimen- 

* = / 0(Am*) s 



sionless isothermal susceptibility Xt = *d( Aft' * f ) A by 



the following equations, 



X. 



qf 



) 2 ~ d (z c ) d r T 



(31) 



As previously mentioned for the critical isochore case, 

XtK) = t r& 1 i while xt(p = i) = K* T {p = 

1) [with xt = (§[s§) A . = 5 ) 2 K rh where k t 
is the dimensionless isothermal compressibility k t — 

There- 



\{%) T ]=^T (With «T = }(g) ). 



fore, the master susceptibility can be also related to the 
dimensionless isothermal compressibility by, 



X, 



(a; 



, 2-d 



(32) 



The master asymptotic singular behavior of X* ( reads as 
follows 



^qf - Z> x 



(T* 



1 



(T* 



(33) 



where 7 = 1.2396935 (author?) [1J|. The master values 
of the leading and confluent amplitudes are Z+ = 0.119 
and Z^ ,+ = 0.555, respectively, where the universal value 



of the confluent amplitude ratio 



0.67919 is given 



in Ref. (author?) [lfj]. The associated asymptotic sin- 
gular behavior of the isothermal compressibility reads as 
follows 



,exp 



(At*) =r + (AT*)~ 7 1 + a+ (At 



(34) 



The term to term comparison of (master) Eq. ([33)1 and 
(physical) Eq. l(34|) . leads to the following amplitude es- 
timations 



(a; 



-7 -7+ 



Z}+ 



(Yc) 



A 



(35) 
(36) 
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Figure 1: (Color online) Singular behavior of the isothermal compressibility of the one-component fluids, (a) kt as a function 
of T — T c > (log-log scale), along the critical isochore for Xe, Kr, 3 He, SFs, CO2, and C2H4 (see inserted Table for fluid 
color indexation); (b) Log-Log plot of tC T (At*); black full curve: mean-field behavior of equation 4 vdw = | (Ar*) _7MF with 

7mf = 1. (c) Matched master behavior (log-log scale) of the renormalized susceptibility X* f — (A* e ) 2 d Z c k* t [see Eqs. I0TJ], 
as a function of the renormalized thermal field T* [see Eqs. i|22]l ]: black full curve: Eq. (J84J) ; red dashed curve: tangent of Eq. 
(jB 1 2 h at the point M (see text and Appendix B); (full) arrow (label PAD): master extension of the Ising like preasymptotic 
domain of Eq. (| 100p : (dashed) arrow (label EAD): effective extension of the extended asymptotic domain of Eq. I|llip 
corresponding to £* f = (A* e ) 1 > 3 (see Ref. (author?) [13]); grey area: master correlation length range 10.5 < £* ( < 0.73 
(thermal field range 1.9 x 10 -2 < T* < 1) discussed in Appendix B . 



with 

Z u+ a + 

= = 0.67919 (37) 

The expected asymptotic collapse of the fluid proper- 
ties on a single curve due to the scale dilatation method 
is illustrated in Fig. Q] (log-log scale). The raw data are 
reported in Fig. to easily distinguish between singu- 
lar behavior of kt (expressed in Pa -1 ) as a function of 
T — T c (expressed in K), for each one-component fluid 
(see the fluid color indexation inserted in Fig. \Vp). Fig- 
ure QJd illustrates the differences between the correspond- 
ing dimensionless behaviors k t (At*) which confirm the 
failure of results provided by the two-parameter corre- 
sponding state principle. This figure also shows the fail- 



ure of mean-field like behavior predicted from the van der 
Waals (vdW) equation of state which is here represented 
by the black full curve of equation K Tvdw (Ar*) 7vdw = 

FydW = g, with 7 v dw = 7mf = 1- On the other hand, 
Fig. \Vp demonstrates the collapse of X* f (T* ) on a mas- 
ter curve where the scatter corresponds to the estimated 
KT-precision (5-10%) for each fluid. We underline the 
combination of the "scaling" and "extensive" roles of the 
characteristic factor Z c in the renormalization [see Eqs. 
(|3"TTl and ([32| ] of the ordinate axis of Fig. [j"b ( compare 
for example with Fig. 3 of Ref. (author?) [2fj| or with 
Fig. 2 of Ref. (author?) [22]). The complementary ma- 
terials for complete analysis of this Fig. Hfc will be given 
below and in Appendix B. 

Therefore, adding the correlation length results given 
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in Fig. lc of Ref. (author?) [33] to the present isother- 
mal susceptibility results, we close the asymptotic mas- 
ter behavior generated by the scale dilatation method, in 
conformity with the two-scale-factor universality of the 
Ising-like systems. 

To summarize the main interest of this Ising-like mas- 
ter description of the {l/}-subclass (associated to the 
selected set {A; v\ 7 } of three independent universal ex- 
ponents (author?) I 1 1 ) . we introduce 

i) the physical amplitude set 



s A = {a+ ; e;r+} 



(38) 



which characterizes the physical Ising-like universal fea- 
tures of each selected pure fluid having the critical set 

ii) the corresponding scale factor set 



Ssf — {Ye, Z c ;A* e } 



(39) 



which characterizes the dimensional universal features of 
the critical interaction cell of each selected pure fluid hav- 
ing (/3 C ) _1 and a c as energy and length units, respec- 
tively, and 

iii) the master amplitude set, 



3. MASTER MODIFICATIONS OF THE MEAN 
CROSSOVER FUNCTIONS 

3.1. System-dependent parameters of the mean 
crossover functions. 

For the {$3 (l)}-class, the mean crossover functions 
F P (t, h — 0) describing the crossover behavior of the the- 
oretical properties P t h (t) as a function of the renormal- 
ized temperature- like field t, for zero value of the external 
ordering (magnetic-like) field h, are given in detail in I. 
All the theoretical functions Fp (t) have the same func- 
tional form whatever P t h ; and, as noted in I, a closed 
presentation of their universal features, only needs to use 
for example the mean crossover functions Fg (t) = j^j^ 

for the inverse correlation length, and F x (t) = Xt *( t ) ^ or 
the inverse susceptibility, at h = 0, in the homogeneous 
phase T > T C [T (T c ) is the temperature (critical temper- 
ature)]. These two theoretical functions read as follows: 



[Ah (*)]■ 



z; 



x+.t D ^ 



(42) 



[X th (tT 1 = Z+ffl (l + X+^Vy*'' (43) 





f 4- + = 


0.555 






0.570481 




{ k= 


0.119 



(40) 



which characterizes the master Ising-like universal fea- 
tures of the {l/}-subclass. Three independent relations, 
i.e., [Eqs. l(28|) . (|35| . and (|36| ]. connecting these three 
previous sets, can be written in the following condensed 
functional form 



S 



{1/} _ 



{S a F(S S f)} 



(41) 



where the function J 7 (Ssf) takes an universal scaling 
form of the two (fluid-dependent) scale factors Y c and 
Z c . Accordingly, any physical amplitude of any one- 
component fluid can be estimated from the equations 
given in Table [J satisfying to the two-scale-factor uni- 
versality of the {$3 (l)}-class (where xenon acts as a 
standard critical fluid to estimate three characteristic 
master amplitudes labeled with an asterisk, see Refs. 
(author?) [3(| [13, 5H]). However, the effective exten- 
sion range where the master behavior is observed, as an 
explicit criteria which defines the preasymptotic range 
where the two-term Wegner-like expansion is valid, re- 
main not easy to estimate precisely only using the scale 
dilatation method. These two problems can be solved us- 
ing a master modification of the mean crossover functions 
(author?) [1.3;] obtained from the massive renormaliza- 
tion scheme, as shown in the next section. 



D (t) is a universal mean crossover function for the con- 
fluent exponents A and Amf which reads 



D (t*) = 



such that D (t = j^yt 



A MF S 2 Vt + A 



s 2 Vt + i 



(44) 



All the universal 



exponents v, 7, A, Amf, and the parameters Z^ , X^, 
Y^, Z+, X+j, and S2, are given in I. 

The temperature-like field t is analytically related to 
the physical dimensionless temperature distance 



At* = 



T-Tr 



by the following linear approximation 
t = At* 



(45) 



(46) 



which introduces $ as an adjustable (system-dependent) 
parameter. Here i? is a scale factor for the tempera- 
ture field. Correlatively, it is important to note that 
the definition of At* [see Eq. jH| ], introduces the crit- 
ical temperature T c as a system-dependent parameter. 
Then the relation between the dimensionless thermody- 
namic free energies of the <I> 4 -model and the physical 
(one-component fluid) system, only involves the energy 
unit {PcV 1 = k B T c . 

Similarly, the ordering-like field h is analytically re- 
lated to the corresponding physical dimensionless vari- 
ables Ajj, (or Ah*) by the following linear approximations 
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(including quantum effects) 

h = tp p 
or h = ip 



(A* qe ) 2 Ah* 



(47) 



which introduce ip p ( or ip) as an adjustable (system- 
dependent) parameter. ip p , respectively ip = (Z c ) 1 ip p , 
is a scale factor for the ordering field A/t, respectively 
Ah* = Z c AJx. 

Accordingly, the dimensional analysis of each term of 
the dimensionless hamiltonian of the $ 4 -model leads to 
the introduction of a finite arbitrary wave- vector Ao, so- 
called the cutoff parameter, which is related to the finite 
short range of the microscopic interaction (see for exam- 
ple, Ref. (author?) [jj). Since the value of the cut- 
off parameter of a selected physical system is generally 
unknown, a convenient method at d = 3 consists in re- 
placing Ao by go (author?) which is the critical 
coupling constant having the correct wavenumber dimen- 
sion (see our introductive part). This system-dependent 
wavenumber go provides the practical "adjustable" link 
between the theoretical dimensionless correlation length 
(•^th) and the measured physical correlation length (£exp) 
of the system at d = 3, through the fitting equation : 



(A* e ) ^exp (At*) ^ (go)- 1 £ th (t) 



(48) 



In Eq. (j48j) . (<7o)~ appears as a metric prefactor for the 
theoretical correlation length function. From Eqs. (@6|), 
(|47)l . and (|48"ll . the asymptotical non-universal nature of 
each physical system is then characterized by the scale 

factor set (.go) 1 ; ipp ( or "^)| (with implicit knowl- 
edge of T c and A* e ). However, for the present fluid study 
where the thermodynamic length unit is already fixed 
by Eq. ©, the above fitting Eq. (|48|) introduces one 
supplementary dimensionless number defined such as: 



goa c 



(49) 



where the notation iS 1 ^ anticipates a master nature of 
this product which we will demonstrate below [see Eq. 
(|93l) ]. More generally, in order to maintain unicity of the 
length unit in the dimensionless description of the sin- 
gular behavior, any theroretical density property (which 
implicitely refers to the length scale unit (go) 1 ) needs to 

introduce the proportionality factor (L^ 1 ^) to the cor- 
responding dimensionless physical density which refers to 
the length scale unit a c . As a direct consequence of the 
fitting Eq. (|48|) for the correlation length, the order pa- 
rameter density m must be analytically related to the 
corresponding physical dimensionless variables Ap (or 
Am*) by the following linear approximation (including 
quantum effects) 



m =(L{W)- d (^)- X [A* e Ap] 
•m = (LW>) d ip- 1 [A* e Am*] 



(50) 



For simplification of the following presentation, we only 
use ip p related to the practical dimensionless form of the 
variables (see above § 2.3). 

Finally, adding the knowledge of the energy unit and 
the length unit for each pure fluid to the theoretical re- 
sults obtained from the massive renormalization scheme, 
the dimensionless singular behaviors of the fluid proper- 
ties are now characterized by the set 



3>SF 



(51) 



made of three dimensionless scale factors (admitting that 
(Pc)- 1 , oi c , and A* e are known). Therefore, it is easy to 
analytically define these three dimensionless parameters 
which characterize each Ising-like fluid, thanks to the ex- 
act values of the mean crossover functions within this 
preasymptotic domain. 



3.2. Three scale-factor characterization within the 
Ising-like preasymptotic domain. 

As already mentionned in the introduction and dis- 
cussed in a detailed manner in I, this asymptotic charac- 
terization is valid within the Ising-like preasymptotic do- 
main where the complete crossover functions of Eqs. f42|) 
and (|43|) can be approximated by the following restricted 
(two-term) Wegner-like expansions (author?) (TH]: 



CPAD.th 



(t) = (z; 



r 



l+Z^' + t A 



PAD,th 



(t) = {z+y 1 t-^[i + z l x +t^ 



(52) 



(53) 



In Eqs. (HJI and {531), z \ + [ see below Eq. {54])], is the 
amplitude of the first confluent correction to scaling for 
the correlation length, which is related to the one for the 
susceptibility 2^'+ [see below Eq. (|55| ]. by the universal 

ratio ^ = 0.67919 (author?) [ill, with: 



i=l 



E x + y+ 



(54) 



(55) 



The theoretical field extension t < Cp 1 ^ of the Ising-like 
preasymptotic domain where the restricted Eqs. {52|) and 
([53]) are valid is defined in I, such as 



/.Mug 
'-pad 



10" 



(S2Y 



1.9 10" 6 



(56) 



Now considering all the theoretical functions estimated 
for all the singular properties of the Ising-like systems 
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(see I), we can note that the universal features in the 
Ising-like preasymptotic domain are characterized by the 
set 



,{MR} 



!}>+ 
-i 



8.56347 

0.471474 

0.269571 



(57) 



of three theoretical amplitudes associated to the set 
{A; 2/; 7} of three universal exponents selected as inde- 
pendent. Accordingly, the restricted forms of two inde- 
pendent fitting equations are 



one-parameter crossover modeling. That infers Ising- 
like equivalence between all estimated crossover functions 
only using three model-dependent characteristic num- 
bers. This result is shown in Appendix A, considering the 
asymptotic crossover infered by the minimal-subtraction 
renormalization scheme (author?) (ill. |24| and the phe- 
nomenological approach given by a parametric model of 
the equation of state (author?) [30j. 

Obviously, from Eqs. (|46l> and l|56p . it is easy to define 
the extension range 



^ ''-PAD 



/-Ising 
■'-PAD 



1.9 x 10~ 



(64) 



z+ 



(tfAr* 



- Z*' + (iJAr* 



(58) 



l + Z^+( l ?Ar*) A 



(59) 

where £ exp and K^.exp are gi ven by the restricted Wegner- 
like expansions of Eqs. <[27|l and ([34]) . respectively. That 
provides the following hierarchical relations 



A* <&- 



(60) 



(61) 



(62) 



with = = 0.67919 (author?) hfl, llj. We un- 

derline the fact that Eq. (|60|) (or equivalently equation 
a~£ = Z^ + d A in the correlation length case), is to be 
first validated (to confer unequivocal Ising-like equiva- 
lence between the first (system-dependent) scale factor 
■d and at). Then Eq. l[6Tj) fixes the asymptotic ampli- 
tude of the dimensionless correlation length and gener- 
ates a single scale factor attached to the selected (phys- 
ical) length unit, which is then mandatory common to 
the thermodynamic and correlations functions. Finally, 
the validation of Eq. f62|) provides unequivocal Ising- 
like equivalence between the second (system-dependent) 
scale factor i\> p and T + (accounting for "critical" and "ex- 
tensive" nature of the susceptibility). 

Equations l|60p to (|62|) satisfy the following condensed 
functional form 



^{MR} 
5 A 



{S A F(§ SF )} 



(/3c 



(63) 



where the function F takes an universal scaling form 
of the dimensionless asymptotic scale factors d and ip p . 
The universal character of Eq. ([63)1 occurs for any 



of the Ising-like preasymptotic domain of the selected 
fluid (labeled with superscript /). Therefore, for each 
one-component fluid for which (or equivalently one con- 
fluent amplitude among a+ or at) is an unknown pa- 
rameter, the remaining question of concern is: How to 
define the validity range At* < £p AD where the theo- 
retical Ising-like characterization by three scale factors 
can replace the experimental characterization by three 
asymptotic amplitudes? 



3.3. Three free-parameter characterization beyond 
the Ising-like preasymptotic domain 

As noted in Ref. (author?) [Ifj, in the absence of in- 
formation concerning the true extension of the Ising-like 
behavior for a real system belonging to the 3D Ising-like 
universality class, the introduction of the scale factors i9, 
tpp, and the wavelength unit go throughout Eqs. f46|) 
to l|48|) cannot be easily controlled. Alternatively, it was 
proposed to introduce three adjustable dimensionless pa- 
rameters Lq £ , Xq £, and using the following fitting 
equations: 



4*xp(Ar*) 



nf =1 



L+(At*) u 



(65) 



(At 



= (x3 ijC )- i z+(Ar*y 



n 



A" 



(t) 



D(t) 



with 



t = ScAt* 



(66) 



(67) 



Lq £ and Xq £ are two adjustable metric prefactors (with 
same value above and below T c ). $£ is a global crossover 
parameter in a sense where it is attached to an unknown 
effective parameter CJ which measures the extent of fit- 
ting agreement involving an undefined number of terms 
in the Wegner-like expansion (see I for details). The de- 
termination of $£ is then equivalent to the determination 
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of & . However, within the Ising-like preasymptotic do- 
main, the restricted forms of the fitting Eqs. ((65)1 and 
(1661) are 



acts as a simple factor of proportionality for the corre- 
sponding leading amplitude (£ + or T + ). The following 
set 



^iE = L* x (z+) 1 (At*)~ v [l + l\ + (tf £ Ar*) A 



exp 



S ijC (Z+)- 1 (At*H [l + Zj?+ (^Ar*) 



(68) 



(69) 



Therefore, the physical leading amplitudes can be calcu- 
lated using the (independent) equations: 



(TO) 



{i? £ ;L5 )£ ;X5 )£ } 



(77) 



is equivalent to the previous set S^f of Eq. lj5lj) . ex- 
cept that the subscript 1C, £ recalls for a single crossover- 
parameter obtained over an extended temperature range 
CJ > £pAD' beyond the Ising-like preasymptotic domain. 
The following condensed functional form 



= {S A ¥ C (Sic, £ )} 



(/3 c )- I ,a c ,A« 



(78) 



i.e., without explicit reference to dc (however the sub- 
script C recalls for the implicit §c dependence due to 
the fitting in the temperature range At* < C* , with 
& > £p AD ). Noticeable distinction occurs for the con- 
fluent corrections to scaling since the first confluent am- 
plitudes are only ^-dependent and can be calculated 
using the equations: 



4 = WcY 



(tf£) A Z 



1,+ 

X 



(72) 



(73) 



0.67919 (au- 



interrelated by the universal ratio -f- 
thor?) [13]. 

For better understanding of the scaling nature of the 
analytical transformations of the physical variables [such 
as Eqs. ([46]) or (J67J)] , we select Eq. ((T3j) as the indepen- 
dent equation for the critical crossover characterization. 
We must then rewrite the above Eqs. |70|) to {73)l in the 
following hierarchical forms 



L 



( x S,r) 



z 



universal est. (74) 



universal est. (75) 



universal est. (76) 



where the l.h.s. of the above equations contain all the 
system-dependent information, first for Ising-like criti- 
cal crossover, then for asymptotic behavior of correlation 
functions, and finally for asymptotic behavior of thermo- 
dynamic functions. Moreover, this information is given 
in a dual form, i.e., as a product between a "physical" 
amplitude (a+, £ + , or T + ) and either a "crossover" factor 
(fie), which acts as a scale factor for the confluent correc- 
tion contribution, or a "pre"-factor (LjJ c or Xq c 



can be used in a equivalent scaling manner to Eq. (|63j) 
(71) when the crossover parameter §c is unique within the & 
range. 

To our knowledge, the unicity of the crossover param- 
eter along the critical isochore of a one-component fluid 
has never been directly evidenced from the singular be- 
havior of the correlation length or any other thermody- 
namic property. However, from simultaneous fitting anal- 
ysis of several singular properties of xenon and helium 3, 
an indirect probe of a single value for one adjustable pa- 
rameter related to the scale factor d was obtained, using 
the crossover functions estimated in the massive renor- 
malization scheme (author?) [l2l l38| and the minimal- 
subtraction renormalization scheme (author?) [H, HH- 
But these results were never used to accurately analyze 
the expected equivalence between Eqs. ((63)) and (f78|. 
and then to estimate the other two scale factors h^ 1 ^ and 
ip P , which is the only correct way to verify the asymptotic 
condition -§ = ■&£ within the Ising-like preasymptotic do- 
main (author?) (47) ]. 

An analytic determination of made beyond the 
Ising-like preasymptotic domain without use of any ad- 
justable parameter, is under investigation for the case of 
the isothermal compressibility of xenon (author?) [4Sl |. 
The main objective is to carefully correlate the local value 
of this crossover parameter with the local value of the cor- 
relation length before to validate its uniqueness by identi- 
fication with the asymptotic scale factor d, calculated by 
using Eq. lf46|l . However, such a challenging demonstra- 
tion of 1? = in the temperature range At* < £g AD , 
i.e., within the so-called Ising-like extended asymptotic 
domain (EAD) in the following, as a formulation of the 
three-parameter characterization of xenon selected as a 
standard one-component fluid, remain two preliminary 
attempts to test the equivalence between Eqs. f63|) and 
((78)1 . That needs to be examinated using a more general 
approach, as the one proposed below, where we will intro- 
duce three master constants which relate unequivocally 
dimensionless lengths and relevant fields of both (theoret- 
ical and master) descriptions, to identify the theoretical 
crossover of the {$3 (l)}-class with the master crossover 
which of the {l/}-subclass. 
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3.4. Identification of the theoretical and master 
asymptotic scaling within the Ising-like 
preasymptotic domain 



Now, while reconsidering our previous analysis of the 
relations between physical and master properties, we 
must rewrite Eqs. l(28|) . (|35| . and (|36l) . in the following 
hierarchical forms 



Y c (a+y 



master est. 



(79) 



3.5. Master modification of the theoretical 
crossover for the isothermal susceptibility 

We start with the following modification of Eq. lf66|) 

N 



=4 1/>a £* 7 II I 1 +•*?*« 



and the following modification of Eq. p7l 

t = e {1/} \t*\ 



,£>*(*) 



(84) 



(85) 



~{YcT 



(A 



qej 



So 



= Z^ = master est. (80) 



by introducing the prefactor Z\ ' and the scale factor 
Q^ 1 ?} as master (i.e. unique) parameters for the {1/}- 
subclass. We note that Q^ 1 ^ , characteristic of the (criti- 
cal) isochoric line (with same value above and below T c ), 
reads as follows 



Z c (Y c y (A* qe ) 2 - d T 



Z + 



master est. 



(81) 



Q{1/} 



rl± 



1,± 



(86) 



Comparison of Eqs. (|74|) to J76J) with Eqs. ((T9j) to (|81j) . 
shows that their r.h.s. differences only concern the re- 
spective numerical values of the characteristic master set 



■f 1 f \ 

S\ of Eq. l(40|) , and universal set 



l { f R} ofEq. (EH. 
For their l.h.s. comparison, neglecting the quantum cor- 
rections in a first approach (i.e. fixing A* e = 1), the term 
to term identification between measurable amplitudes 
underlines the analogy between the explicit parameter set 
{Y C ;Z C }, related to the master description, and the im- 
plicit one {"!?; ip p }, related to the massive renormalization 
description. We can then note that the {l/}-master for- 
mulation compares to the $3 (l)-universal formulation, 
only if we have correctly accounted for the asymptotic 
scaling nature of each dimensionless number needed by 
the massive renormalization scheme. In order to reveal 
such a scaling nature, it is essential to note that the scale 
dilatation method replaces the renormalized fields (such 
as t, h, to, etc.) needed to observe the "universal" be- 
havior of the $| (l)-universality class, by the {l/}-fields 
(such as, T* , H^f, A4* f , etc.) needed to observe the 
"master" behavior of {l/}-subclass. The common physi- 
cal variables are At*, A/2, and A/5. Therefore, it remains 
to give explicit forms for the following exchanges between 
the theoretical variables and the {l/}-subclass variables 



whatever the selected one-component fluid is. By virtue 
of the universal feature of confluent amplitude ratios (see 
Table I), the numerical value 



6 {1/} = 4.28810" 



(87) 



is the same whatever the property and the phase domain. 
However, we also note that Q^ 1 ^ contributes to the lead- 
ing term. Thus, in addition to Eq. (|84|) . we define Z 



{1/} 



such that 



4 1/} 



The numerical value, 



Z[ lf} = 1950.70 



(89) 



is the same in the homogeneous phase and in the non 
homogeneous phase. The curve labeled MR in Figure Q] 
was obtained from Eqs. |84|) and (|85)l using the numerical 

values of 9 {1/} and z[ lf} given by Eqs. p? ) and ([89 ]l . 
respectively. 

We recall that our previous analysis (author?) [39j of 
the correlation length has introduced a similar prefactor 



{1/} 



through the following modification of Eq. (J65 



t -> T* 



H 



qf 



M 



qf 



(82) 
(83) 



(see Ref. (author?) [33] for the correlation length case). 
The next subsection is dedicated to the isothermal sus- 
ceptibility case (which then closes the description of the 
{l/}-subclass along the critical isochore in conformity 
with the universal features estimated for the Ising-like 
universality class). 



with 



%l lf} F e (t) 



zfzf 



which has the same numerical value 
Z\ lf} = 25.6936 



(90) 



(91) 



(92) 
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(a) F P 

i { A MR} \\ m 



00 



°2P1S 



■{!/} 
SF 
{1/} 
A 



(c) 



SsF 

S A 



'P. 



V 



"qf 



^qf 



Hoi) 

l|40|l 



p* 

exp 



dzzl 



Kj 1 



dm 

(f39T> 

(EH) 



W 7 " 



0.471474; 



0.269571 



= 25.6936; l{ in = 1950.7 



7 wy 



h^f} = 25.6936 



(Li 1 f}) d = 16961.9 



1950.7 



6 C _ .,,±/ (lfl\2- 



«z±(e{i/l) 
z 



^M—j, = (LlWf = 2.94878 



{lW} = 25.6936; *W> = 1.73847 10" 4 } 



i 



27 = 0.570481; = 0.119 



i 



{h* ) c ; Xg £ } 



= 1.20999A* e (y c )-' 

■ Z c (Y c p 
: 0.44144 



( L f 1 /}) d (vi/{i/}) 2 (e{i/})- 



= (A; e )- d (y c )2-^| =(A; e )- d (Y c )2-(L{i/}) d (0{i/}) 5 
= 0.564481 (A; e )" d (y c ) 2 -° ! 

= ( y ^ g i |y = golf ; (L {i/}) d vi>{i/} (ed/})^ 



A* e (Z c )5 

0.499185- 



= 25.6936; V P = (Z c )~i -P* 1 /}} 



{ 



25.6936; A! 



} 



7l,+ 



^ 

= 8.56347 



2 



i.± 



0{i/} _ 4,288 x 10~ 3 



9 {i/} = 
{i/} = 

Q{1/} = 
©{I/} 



= 4.288 10 -3 
= 0.555 



■a L = # = y c eW> 



= y c ew 



Table II: Three parameter characterization (colum 3: leading amplitudes or prefactors; colum 4: scale factor or crossover 
parameter; see text): (a) for the mean crossover functions Fp (t) defined in I; (b) for the master crossover functions "P* f (T*) [see 

Eq. ((109 }) ] ; lines 4, 9 and 10: independent parameters; lines 5 to 8): related parameters; The two relations (iS^^ d = Z^ 1 ^ 
are in conformity with the two-scale-factor universality; (c) similar to (b) for the physical crossover 

functions P e * xp (At*) [see Eq. I|110p ]: The two relations ( ) d Ljj £ = 1 and (^o,c) d 7 X °' C yi = ^ are m conformity with 

the two-scale-factor universality; All the values of c , $c = $ and ipp can be estimated from Q" lln = {(/3 C ) _1 ,a c ,Y c , Z c } and 
A* e of the selected one-component fluid. 



and 



(4 m r 



7TTTT 



= Z 



for the homogeneous and non homogeneous domains. Of 
course, we retrieve here the previous Eq. lj49f 



Z 



Vfluid 



(93) 



which now is valid whatever the fluid under considera- 
tion. The set of master (two pre- + one scale) factors 



S, 



{if} 

2P1S 



6^ 1/ > = 4.288 x 1(T 3 
{ij 



Z { c lf} = 25.6936 



(94) 



Z 



1950.70 



closes the universal behavior of the {l/}-subclass, as 
shown by the results reported in Table HT1 for all the prop- 
erties calculated along the critical isochore (for notations 
see below and Refs. (author?) (HI, IH, 51). Equation 
(|49| [or Eq. (|93| ] appears then as the basic hypothesis 
which defines the critical length unicity (author?) [ItJ 
between correlation functions and thermodynamic func- 
tions of the one component fluid subclass. HJ 1 -^ takes 
an equivalent nature to the length reference used in the 
renormalization scheme applied to the $3 (l)-class, what- 
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ever the selected physical system. 

The major interest of Eqs. lf88| and lf9lj) is that 
they introduce the needed "cross-relation" between pure 
asymptotic scaling description and first confluent correc- 
tion to scaling, in order to obtain only two independent 
leading amplitudes within the Ising-like preasymptotic 
domain. Such a cross-relation occurs if the non-universal 
scale factor associated with the irrelevant-field which in- 
duces the correction-to-scaling term of lowest relative or- 
der (Ar*) A in a Wegner-like expansion, is the same as 
the non-universal scale factor associated with the rele- 
vant (thermal) field which gives the leading scaling term 
(Ar*p p . 

In that universal description of the confluent correc- 
tions to scaling, each crossover function includes the 
(two-term) master behavior expected using the scale di- 
latation method. By comparing the leading terms on 
each member of Eqs. ([84|. (|33)l . and 1(69]). we obtain the 
relations 



r ± = 



(z^>z±) V^=x* £ ( 



(95) 



where the fluid-dependent metric prefactor Xq c of Eq. 
(l69l) now reads as follows 



a c = z±z± (Y c y 



(96) 



In Eq. ([96|) , the critical contribution of the scale factor Y c 
is explicit. The remaining adjustable crossover parameter 
i?£ of Eq. f67j) is characteristic of the Ising-like extended 
asymptotic domain At* < £g AD where the theoretical 
crossover functions and experimental data agree. Within 
the Ising-like preasymptotic domain [see Eq. lj56|) ] where 
the two-term Wegner-like expansions are expected to be 
valid, the comparison of the first confluent amplitudes for 
master and theoretical descriptions, enables one to write 
•dc as follows 



(= r c ew> 



(97) 



3.6. Closed master modification of the mean 
crossover functions and master extension -Cj^^, of 
the extended asymptotic domain 



Obviously, the equivalent approach at exact critical- 
ly and along the critical isotherm occurs in virtue of 
the two-scale-factor universality which implies a second 
unequivocal relation between ip p and Z c . However, we 
can anticipate such a result only from the thermody- 
namic definitions of the susceptibilities xth = (^jt) t 



= ( 



ax- 



introducing the scale factor ty^H 



through the following linearized equations 



h = W>H* f = (A, 



(*{i 



M 



qfl 



A* e \M* 



(101) 



(102) 



where vpf 1 /} is a master (i.e. unique) parameter char- 
acteristic of the (critical) isothermal line for the {1/}- 
subclass (vfff 1 /} has the same value whatever the sign 
of the order parameter). From comparison between ei- 
ther Eqs. m|, (El, 021 and (EH) or Eqs. (2D, flU, 
(l49l) and l|102j) . it is immediate to show that Xth = 



(L* 1 ^) d (*W>) 1 X* t and, correlatively, to obtain the 
following expected relation 



tP p = (Z c ) Mi/1 



(103) 



The unequivocal link between the scale factors needed, 
either by the theoretical description, or by the master 
description, is given by Eqs. {93]), J97|) and (|103j) . There- 
fore, the leading theoretical and master amplitudes of the 
susceptibility and the order parameter are related by the 
equations : 



with 



(104) 



4'± 



Z 



7-* (— Qi 1 /} vj* 



(98) 



1m 



(105) 



As considered from basic input of the scale dilatation 
method, Eq. (|98|) agrees with the scale dilatation of the 
temperature field 



T* = Y c I A? 



(99) 



— i 1 f i 

Note that the extension T* < £p A p of the Ising-like 
preasymptotic domain of the {l/}-subclass can then be 
immediately obtained from Eq. (|56j) . with 



C 



{1/} 

PAD 



/rising 
^PAD 



4.43 x 10" 



(see for example the full arrow labeled "PAD" in Fig. . 



while the leading theoretical and master amplitudes of 
the correlation length and the heat capacity are related 
by the equations : 



zfzf 



Z ± 



1/} (onny 



c 



\ d / , , \ 2 — a 



(106) 



(107) 



(100) where the master prefactors Z^ 1 ^ an ^ are f° r the 



7{1/} 



heat capacity case and the order parameter case, respec- 
tively [see below, Eq. (j 109[) ] . Finally, the characteristic 
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MR master fluid 

Sa R =^a ] ^{S > sf') S [ a' = S a F(S sf ) 
s MB - w c(if) - w .<? . 



ST 



Figure 2: Schematic links between three amplitude characterization §^ Mi ^ of Eq. I|57p . 5^^' of Eq. l14()[l and <Sa of Eq. 
(l38l> . of the theoretical, master and physical singular behaviors, respectively, for a fluid of critical parameters given by Q™ m = 
{(/3c) -1 , c*c, y c , Z c } of Eq. ((3j) and A* e , and which belongs to the one-component fluid subclass. 



set 

feW* = 4.288 x 1(T 3 ) 

LW> = 25.6936 \ (108) 

= 1.73847 x 10~ 4 J 

is Ising-like equivalent to the one of Eq. f94|) and 
closes the modifications of the theoretical functions of 
the $3 (l)-class in order to provide accurate description 
of the master singular behavior of the {l/}-subclass. 
Accordingly, each modified function reads as follows 

V* f (T*)=Z { p lf} F P (t) (109) 

with t = e^T* and F P (t) defined in I. All the master 
prefactors Z p can then be calculated using the rela- 
tions given in part (a), column 3, of Table HT1 Within 
the Ising-like preasymptotic domain, Eq. (|109[) can be 
approximated by Eq. Q30p . 

Alternatively but equivalently, each physical property 
can also be fitted by the following modified function 

N y± 

P e ; p (|Ar*|) = P^Zp- |Ar*r P LI ( l + X tp tD{t) ) 

(110) 

with t = | At* | = QWlYc \At* \ and where the function 
D it) [see Eq. Q44p ] and the universal quantities Zp, ep, 
Xf p , Y^p, are given in I. All the physical prefactors P*, c 
can also be calculated using the equations given in part 
(b), column 3, of Table ITU where the physical prefactors 
Cox an d M ,£ are for the heat capacity case and the 
order parameter case, respectively [see Eq. QUO)) ]. 

As a summarizing remark related to the schematic Fig. 
[21 the theoretical amplitude set S A MR ^ of Eq. (|57]l . the 

-f 1 f \ 

master amplitude set S A of Eq. |40|) , and the physical 
amplitude set Sa of Eq. (|38|) , are unequivocally related 
only using Y c and Z c (or t9 [see Eq. l(97j) ] and ?/> p [see Eq. 
Q103P ]) as entry parameters (assuming that ((3 C )~ , etc, 
and A* e are known). 

In addition, we can also account for the results of pre- 
vious analyses of different singular properties for several 
one-component fluids where each master singular behav- 
ior is well-fitted by the corresponding crossover functions 



in the extended asymptotic domain which corresponds 
to t* i{ > 3 — 4 (see for example the dashed arrow labeled 
"EAD" in Fig. for the susceptibility case). Indeed, 
the effective extension where this modified the- 

oretical description seems to be valid, corresponds to the 
temperature-like range such as 



T* < C^ f} ~ 0.07 - 0.1 (111) 



Equations QlOOp and Qllip are of crucial importance for 
experimentalists interested on liquid-gas critical point 
phenomena since they are the "master" (experimental) 
answer to the unsolved theoretical question: How large 
is the range in which the asymptotic universal features 
are valid in pure fluids? Moreover, when Q™ ln = 

|(/3 C ) 1 ,a c ,Y Cl Zc^ and A* e are known, we note that 

each modified crossover function of Eq. <| 109[) can act 
beyond the Ising-like preasymptotic domain, i. e., within 
the two-decade range 10~ 2 < T* < 1 corresponding to 
the grey areas of Fig. [1] to confirm that the critical Ising- 
like anomalies characterized by a limited numbers of crit- 
ical parameters would dominate in a large range around 
the liquid-gas critical point. Such a modified theoretical 
analysis of the available fluid data at finite temperature 
distance appears then similar to the one initially pro- 
posed to provide the first test of the scaling hypothesis 
for the one-component fluids by using effective univer- 
sal equations of state with only two adjustable dimen- 
sionless parameters. As a typical example, we analyze 
the isothermal susceptibility for twelve different fluids in 
the Appendix B, using the well-known linear model of a 
parametric equation of state (eos) (author?) with 
7eos = 1-19 (and /3 eos = 0.355 to close "thermodynamics" 
scaling laws). Furthermore, Eqs. QlOOp and Qllip of- 
fer explicit Ising-like criteria to control the development 
of any empirical multiparameter equation of state where 
such a minimal critical parameter set QJ? ln is customarily 
used (see for example Ref. (author?) [5(j] and references 
therein) . 
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4. CONCLUSIONS 

We have shown that the needed information to describe 
the singular behavior of one-components fluids within 
the Ising-like preasymptotic domain was provided by a 
minimum set of four scale factors which characterize the 
thermodynamics inside the volume of the critical inter- 
action cell. We have illustrated the Ising-like scaling na- 
ture of the scale dilatation method able to demonstrate 
the master singular behavior of the one component fluid 
subclass. Using the mean crossover function for suscep- 
tibility in the homogeneous phase, which complements 
a previous study of the correlation length in the ho- 
mogeneous phase, we have demonstrated that the uni- 
versal features predicted by the massive renormalization 
scheme is then accounted for by introducing one common 
crossover parameter and appropriate prefactors, only two 
among the latter being fluid-dependent. Defining three 
master constants able to relate the theoretical fields and 
the master fields, the corresponding master modifications 
of the mean crossover functions were obtained from iden- 
tification to the asymptotical master singular behavior of 
the one-component fluid subclass. The four critical co- 
ordinates which localize the gas-liquid critical point on 
the pressure, volume, temperature phase surface provide 
then the four scale factors needed to calculate the sin- 
gular behavior of any correlation function or thermody- 
namical property, in a well-controlled effective extension 
of the asymptotic critical domain for any one-component 
fluid belonging to this subclass, in agreement with the 
idea first introduced by one of us. In the case where 
quantum effects can be non negligeable, a single supple- 
mentary adjustable parameter seems needed to correctly 
account for them. 
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Table III: Estimated universal values of the confluent expo- 
nent (column 2) and confluent "crossover parameter" (column 
3) of the scaling forms of Eq. l|A2[l for the first confluent cor- 
rection term in the susceptibility case. Column 1: label M 
of the different crossover models (see references given in the 
last column). Column 5: order of magnitude for the ratio of 
the crossover parameters obtained using MSR or CPM fitting, 
from reference to the MR fitting (see text). 



Although crossover phenomenon can be general upon ap- 
proach of the Ising-like critical point, such a modeling, 
in which G is a tunable parameter, is essential to check 
carefully its description with the objective to discuss the 
shape and the extension of the crossover curves (lead- 
ing for example to distinguish a wide variety of Ising-like 
experimental systems, including simple fluids, binary liq- 
uids, micellar solutions, polymer mixtures, etc.). How- 
ever, for the one-component fluid case, our interest can 
be restricted to the crossover temperature scale estimated 
by three crossover modeling selected in Table Hilt i.e., i) 
the massive renormalization scheme (labeled MR) (au- 
thor?) [n], 0] and ii) the minimal subtraction renor- 
malization scheme (labeled MSR) (author?) [Til, El, 
both modeling without tunable G, and iii) the paramet- 
ric model of the equation of state (labeled CPM) (au- 
thor?) f^jl, with tunable G. The universal form (au- 
thor?) [13, l2J, |30j of the first confluent amplitude for 
the susceptibility case, is then given by the equation 



9\,M 



(A, 



> A' 



(A2) 



Appendix A: SCALING EQUIVALENCE FOR A 
ONE-PARAMETER CROSSOVER MODELING 
WITHIN THE PREASYMPTOTIC DOMAIN 



The use of Eq. J60J) in the hierachical Eqs. <[60j) to 
([62"ll . needs that the characteristic scale factor d is the 
first mandatory parameter to be determined, whatever 
the renormalization scheme (at h = 0). For scaling un- 
derstanding, Eq. (|60|) must be expressed in the universal 
form of Eq. f74|) . i.e., such as 



It) 



-An 



(Al) 



Such a theoretical scaling form of Eq. (|Al[) [or Eq. 
(|74|) ] is then provided from any phenomenological model 
which use a single crossover (temperature-like) param- 
eter At* m related to the (system-dependent) Ginzburg 
number G (the subscript M refers to the selected model). 



where 9^~m is an universal constant given in Table 

The differences in the estimates of g^j^ account for dif- 
ferences in several theoretical aspects: the extension of 
the renormalization procedures, the nature of the asymp- 
totic limit of A gr- , the nature of the non universal cor- 
rections, the numerical calculations, etc.. Therefore, we 
cannot expect practical understanding from each value 
given in Table Mil However, in spite of these numerical 
differences, the scaling form of Eqs. ([73|. £51]), and (|A3[) 
provides analytic equivalence between the three models 
since each model exactly accounts for the same Ising- 
like critical crossover using a single crossover parameter, 
especially for temperature dependence of the effective ex- 
ponent (author?) [25[. The crossover temperature scale 
Ar* M takes a small finite value and can then be "compa- 
rable" to 1?, via the "sensor" Ar^ = f [see Eq. (39) in I] 
of the mean crossover functions (see also Ref. (author?) 
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Figure 3: (Color online) Schematic illustrations (for 7 — 7e and 

(aj) A ) of the Ising-like preasymtotic equivalence between 
the dimensionless crossover temperature scale At* m needed 
by the one-parameter crossover model M (dashed red line), 
and the scale factor d needed by the massive renormalization 
(MR) scheme (full black line) [see Eqs. {H}, lfA2|) . and lfA3|l . 
and Table Ull] ; The small difference on the respective A values 
is neglected. 



[48j]). As illustrated by the point to point transforma- 
tions in Fig. [3K and b, and numerical values given in 



column 5, Table Hill At* m is then scaled by t? through 
the "universal" scaling equation 




tfAr. 



universal est 



(A3) 



where Ar* MSH = b* + ^-{l-^) & for the minimal 
subtraction renormalization scheme, and At* cpm = 

C* - 1 ' 

TT^p (1 — u) 3 for crossover parametric model, are the 
so-called effective Ginzburg numbers (see the Refs. (au- 



thor?) [24], [3(J for the notations and definitions of the 
above quantities). 

Correlatively but uniquely when Eqs. (|78j) or (|A2j) 
are valid (i.e., when the Ising-like critical crossover is 
characterized by a single parameter), we must extend 
the scaling analysis to the leading amplitudes, expressing 
again Eqs. (j6lj) and (|62j) in the "universal" form of Eq. 
(1551 . i.e. such as : 



(z+)- x = r+[(LW>) d (v,r 2 ^ 



(A4) 



(A5) 



Obviously, as for the confluent amplitude, we can close 
the asymptotic identification between the three (MR, 
MSR, CPM) modeling, introducing two supplementary 
universal numbers which relate unequivocally the scale 
factors h^ 1 ^ and tp p of the massive renormalization 
scheme, to the equivalent two free parameters of another 
crossover approach (see also Refs. (author?) [H, SH| 
and the § B3 below). 



Appendix B: EFFECTIVE CROSSOVER 
FUNCTION BEYOND THE ISING-LIKE 
PREASYMPTOTIC DOMAIN 

1. Effective exponent and effective amplitude 

According to the above asymptotic analysis of the 
equivalence between crossover modeling, the scale trans- 
formations of the variables which produce the universal 
collapse of the Ising-like crossover curves can be illus- 
trated by using, not only effective exponents (author?) 
[25[, but also effective amplitudes (see also Ref. (au- 
thor?) (Hi). Indeed, from xth (*) of Eq. (03j, the local 
value of the effective (theoretical) exponent 7 e ,th (*) 1S 
defined by the equation 



7e,th (*) 



dLn [xth (*)] 
dLnt 



(Bl) 



The local value of its attached effective (theoretical) am- 
plitude Z+ x (t) is defined by the equation 



K.e (*) = £ 



Xth (t) 



(B2) 



Therefore, j e ,th (t) and Z+ (t) have equivalent "univer- 
sal" features as \th (t). By eliminating t [then simultane- 
ously eliminating the scale factor #c since t = ^Ar*], 
the classical-to-critical crossover is characterized by a sin- 
gle (i.e. universal) function Z+ e (7e,th) over the complete 
range 7mf < 7e,th (*) < 7- This result is here represented 
by the top (black dot-dashed) curve in Fig. SJ Its limit- 
ing Ising-like critical point takes "universal" coordinates 

1 7; (Z+) 1 | (see the top cross in Fig. [4|). 
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Eq.(B2) 

Eq.(B6 

Eq.(B4) 
exp. 




Figure 4: (Color online) Theoretical estimations of the effective mean [Zj e , upper black dot-dashed curve, Eq. ((B2[) ] . master 
[Zx,e> median blue double dot-dashed curve, Eq. ((B6[l ] and physical (xenon) [r+, lower red dashed curve, Eq. (|B4[) ] amplitudes 
as a function of the effective exponent 7 e for the susceptibility case along the critical isochore in the homogeneous domain; 
Double (dotted at 7 e = 7, full at 7 e = 7eos) arrays: point-to-point (plusses at 7 e = 7, open circles at y e = 7eos, open squares 
at 7e = 7mf) transformations between effective functions using Y c and Z c , or, alternatively but equivalently, $ and Xg t £ (each 
relation associated to the transformation at 7 and 7 e constant value is illustrated in an attached rectangular box); Lower black 
dashed bold curve (labeled exp): F^from Giittinger and Cannell's fit for xenon susceptibility (author?) [49l | (see also text and 
Table ITV]) : M-coordinates 7 eos = 1.19 and Z^ e = 0.15374: tangent line at the point M to the theoretical curve of Eq. ll84l in 
Fig. [TJ Others quantities, points, and symbols: see text. 



In a similar way, from the physical function 
K Texp (At*) °f Eq. lj66|) which fits the experimental re- 
sults using -&c [see Eq. I(g7]l] and Xg £ [see Eq. ((96j)] , the 
local (physical) exponent is defined by 



7e 



.(At* 



dLr 



(At*)] 



dLn (At*) 

and its related local (physical) amplitude by 



T+ (At*) 



K 



T,exp 



(At*) 



( AT *)-7=,ex P 



(B3) 



(B4) 



Eliminating At* from Eqs. l|B3|l and jB4} , the corre- 
sponding physical function (7 e , e xp) is represented in 



Fig. [4] by the bottom (red dashed) curve, selecting xenon 
as a typical example (author?) [4a |. Its related Ising- 
like critical point takes the physical coordinates {7; 
as represented by the bottom cross in Fig. [4] (with 
r + (Xe) = 0.0578204). For quantitative comparison in 
this "physical" part of Figure HI we also have represented 
the experimental lower (black dashed) curve for T+ val- 
ues obtained from the Giittinger and Cannell's fit of their 
susceptibility measurements (author?) (bold part of 
the curve) , and from several pVT measurements reported 
in Table |IV] (full points labeled 1 to 4, open circle labeled 

P). 

Finally, considering the master singular behavior 
X* { (T*) of Eq. HI using O^ 1 ^ [see Eq. (J67j)] and z[ lf} 
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[see Eq. l[96|) ]. we can define the local (master) exponent 

by 

dLn \X*JT*)} 

and its related local (master) amplitude by 

#* f (T*) 



(T* 



-7e,I/ 



(B6) 



After T* elimination between Eqs. (|B5|) and (|B6|1 . the 
master function Z+ e {j e ,if) can also be represented by 
the unique median (blue double dot-dashed) curve in Fig. 
IH Its Ising-like critical point takes the master coordi- 
nates {7; Z^}' corresponding to the median cross in Fig. 

H 

Our main interest can then be focused on the point 
to point transformation at constant 7 e between these 
three curves, using only two fluid-dependent parame- 
ters, either ??£ and Xj$ £ for the physical quantities, or 

O^f} and Z^ 1 ^ for the master quantities. We recall 

that when d c (respectively 9 {1/} ) and z\ lf} = L^l 
are known, XJj £ gives unequivocal determination of ip p 

(respectively ^f^- 1 ^). Now, introducing also Y c and Z c , 
the complete set of the relations between the - theo- 
retical, master, and physical - amplitudes are summa- 
rized in Fig. [H Consequently, this figure closes the 
master description of <^qf (7e,i/) establishing unequivo- 
cal link between the three parameter sets {$£;Xq £ }, 

jo^ 1 ^; Z^H, and {Y C ;Z C }, and also contains explicit 

equations of the schematic links given in Fig. [2] for the 
isothermal susceptibility case [with the implicit master 

condition Z^ 1 ^ = = goa c fixing go]. 

Hereafter we discuss the experimental results obtained 
at large distance to the critical point, i.e., beyond the 
Ising-like preasymtotic domain where practical estima- 
tions of 7 e are significantly different from 7 (an anal- 
ysis of the Ising-like preasymptotic domain very close 
to the Ising-like limit je — * 7 will be in consideration 
in Ref. (author?) [48J; see also below § B.3). Es- 
pecially we focus our present attention on the range 
1.215 > 7 e > 71 w 1.12 corresponding to the grey area 
in Fig. [4] (obviously equivalent to the grey area in Fig. 

at). 

We start with the xenon (Xe) case selected as a stan- 
dard one-component fluid. We can then estimate the 
(theoretical, master and physical) crossover functions for 
the correlation length and the isothermal compressibil- 
ity of xenon, using {Y c = 4.91373; Z c = 0.28601} 
and {&c = 0.021069; X* c = 0.214492}, (or 

; g = 29.1473 nm- 1 }), 
0.443526 (for detail, 



{1? = 0.021069; V P = 3.2507 x 10" 



with a c = 0.881508 nm and L 



ox 



see Ref. (author?) [43]). As a basic application, we 
can define the correspondence between theoretical and 
physical temperature range and between theoretical and 



physical correlation length range for description of either 
7 — 7e,th and Z+ e — (Z+) as a function of t and as a 
function of £ t h , or 7 — 7 e ,exp and — T + as a function 
of At* and as a function of £* . Each respective result 
is illustrated by a (black dot-dashed or red doted) curve 
in each part a to d of Fig. 0D Now, the grey areas 
in Fig. [5] correspond, either to the theoretical ranges 
<i<2x 10~ 3 (bottom axis) and 180 > £ th > 18 
(top axis) in parts a and b, or the physical (xenon) 
ranges 5 x 10~ 3 < At* < 10 _1 (bottom axis) and 
10.5 > £* > 0.73 (top axis) in parts c and d. 

The following compares these theoretical predictions to 
the 7e,pVT and ^tpVT values obtained from pVT mea- 
surements (author?) IHll [H2I [53L [Bil . [HH (see also de- 
tails in Ref. ( author ?)Q). We recall that the pVT 
measurements were performed at finite distance to the 
critical point, such that the n^pVT data obtained from 
pVT data can be fitted by an effective power law 



K 



T,pVT 



"7e,pVT 



(B7) 



only valid in a restricted temperature range defined by 
Ar* in < At* < Ar* lax . The measured (exponent and 

amplitude) parameters ^f e ,pVT', ^tpvr} are then associ- 
ated to the temperature range {At^; Arj^ ax } of central 
value (Ar* pyT ) = V Ar mm Ar max ( in log scale) located 
beyond the Ising like preasymptotic domain. Therefore, 
we can represent these results by points of respective 

coordinates \j e , P VT\ ^X p vt\i {le, P vr; (At* vt )} and 
jr^yj,; (At* vt )| in each appropriate binary diagram. 

The four points (labeled 1 to 4) illustrated in Figs. 
31 |U and HJi, correspond to the xenon results reported 
on lines labeled 1 to 4, respectively, of Table IIV1 The 
points labeled 1 and 2 follow the general trend of the the- 
oretical curves. This result confirms that, in spite of a 
large correlated error-bar in the adjustable exponent and 
amplitude parameters, the variations of their respective 
central values agree with a two-parameter description 
within the "Ising-like" side of the crossover domain where 
7 > le.pVT > 1.17 . However, the point labeled 3, and 
more significantly the point labeled 4, show that thepVT 
experimental results are not in agreement with the mean- 
field behavior predicted by the crossover function within 
the "mean-field- like" side where 7.1 > j e ,pVT > 7mf- The 
failure of the classical corresponding state theory is also 
illustrated by the point labeled 5 in Fig. 01 which cor- 
responds to the result obtained from the van der Waals 
equation of state [see the line labeled 5 (vdW)) in Table 

El. 

To translate the £eae^~ mas ^ er value [Eq. I|llip ] in a 
7£- master value which delimits the effective range of the 
extended asymptotic domain in Fig. (H one needs to con- 
sider the upper horizontal axis of Figs. [St and [Hi which 
measures the master correlation length £* = - 4^ [i.e. 
the dimensionless ratio which compares the size of the 
critical fluctuation to the actual range of the microscopic 
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Je,pVT 


1 e,pVT 


<Ar p V T ) 


Ar4 (7e,th) 


1 e,th 








/ A * \ /A h— * A 


7e,th (l) = 7e,pVT (Ar ) 




1 


1.211 ± 0.01 


0.0743 ± 0.015 


2.07 x 10 -3 


2.95 x 10 -3 


0.07263 


2 


1.16665 


0.089 


2.24 X 10 


3.338 X 10 




3 


1.1198 (= 7,1 j 


0.101 


1.21 X 10 


i nno 1A — 1 

1.928 x 10 


0.09960 


4 


1 (= 7mf) 


0.11 


7.1 10" 1 


oo 


0.08507 


5 (vdW) 


1 (= 7vdw) 


i (- r+ ) 




oo 


0.08507 


P (eos) 


1.19(= 7eoS ) 


0.0793 (= r+ os ) 


1.13 x 10~ 2 


1.135 x 10~ 2 


0.08084 



Table IV: Column 1: index of the points in Figs. [HandOS Columns 2 and 3: Effective power law description oipVT measurements 
in xenon (see Ref. (author?) [48] for detail and data sources); Columns 4 to 6: calculated values of the (geometrical) mean 
temperature (At* vt ) = sj Ar£~ Ar,* ax of pVT measurements (column 4), theoretical local temperature Ar t * h (7 e ,th) satisfying 
the condition 7 e .th (t) = "f e , P VT (Ar*) (column 5), and theoretical local amplitude r^~ th for 7 e ,th (t) = Je, P vr (Ar*) (column 6). 




Figure 5: (Color online) a) dot-dashed (black) curve (labeled MR): 7 — 7 e ,th as a function of t, calculated from the theoretical 
crossover function of Eq. lf43l for susceptibility (log- log scale); (green) line Ta: limiting singular behavior [see Eq. (|B18[) ] 
within the Ising-like preasymptotic domain of extension t* < Cp'J^ [vertical green line, see Eq. l[56jl ]: curve Ta of slope A 

crossing the vertical line t — 1 (x): value of the first confluent amplitude of Eq. l[55jl : vertical (blue and pink) lines: t 11 

2 

and fA-coordinates for 7e it h = 71 and D (t&) — Ai, respectively, b) dot-dashed (black) curve (labeled MR): same as a) for 

(Zj e ) 1 — (Zj) 1 as a function of t, calculated from Eq. I143P : (green and blue) lines T e and T z : logarithmic singularity and 
power law approximation, [see Eqs. (|B 14[) and (|B15[) . respectively, and text for detail]; c) and d) dotted (red) curves (labeled 
Xe-MR): corresponding xenon quantities 7 — 7e and r+ — T + as a function of Ar* = [i9 (Xe)] -1 i; points P and 1 to 4 : pVT 
results (see Ref. (author?) [48fl for detail) given in Table [TV] (see also Figure [4}; ( green and blue) lines Pa and P z : xenon 
counterpart of the theoretical (green and blue) lines Ta and T z . x at t = Ar* = 1 : point-to-point transformations between 
first confluent amplitudes and a J (upper arrow) and between leading amplitudes (Zj) and r + (lower arrow). 
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interaction, with A* e (Xe) = 1 in xenon case]. As a mat- 
ter of fact, the value J~^\ B s» 2 x 10~ 2 corresponds to 

the value yJlxe) wnere £* ~ 3. Therefore, the associ- 
ated local value is ~ 1.16 — 1.17. This value descrim- 
inates the "non Ising-like" range j e < Jc (including the 
value 71 = 7+ 2 MF ~ 1.12) where the effective classical- 
to-critical crossover for xenon is no longer accounted for 
by the theoretical crossover function, as shown in Fig. 
H] where it is observed an increasing difference between 
the curves labeled Xe-MR and the dotted curve labeled 
Xe-exp when 7e — > 7mf = !• 

Accounting for an extended (Ising-like) asymptotic do- 
main defined by 7 e < 7 £, we are also able to revisit the 
results previously obtained using an universal scaled form 
of the equation of state with "universal" values of the ex- 
ponents significantly different to the "Ising" ones. As a 
typical example, the xenon results obtained from the re- 
stricted linear model of a parametric equation of state 
with 7eo S =1-19 [see the line labeled P (eos) in Table [TV] 
are in excellent agreement with the theoretical crossover 
function, as illustrated by the point labeled P in Fig. 
|H Moreover, using as a Ar*-coordinate the theoretical 
value Ar* h (7 eos ) = 1-135 x 1CT 2 [Eq. <[B3|) ]. we can show 
that these results are also well accounted for in the the- 
oretical temperature dependence (see the corresponding 
points labeled P' and P" in Figs. and[5]i, respectively). 
Such results confirm that two xenon-parameters involved 
in an universal form of the equation of state can be used 
as Ising-like characteristic factors to be related to the two 
scale factors Y c and Z c , as illustrated in the next Section 
for the case of the linear parametric equation of state. 



2. Master crossover provided by a restricted linear 
model of a parametric equation of state 

In the seventies, first analyses of the two-scale fac- 
tor universality for one-component fluids used effective 
scaled forms of the equation of state (eos) to fit the pVT 
data measured at finite distance to the critical point (for 
detail see Refs. (author?) 0, [jj, [HE [H3, \M, IH 0, 
[6l|). Such a thermodynamic approach of universality 
was based on a limited number of characteristic param- 
eters for each pure fluid, using effective universal values 
for the critical exponents. We limit the present purpose 
to the well-known restricted linear model of the paramet- 
ric equation of state (author?) [61], with application to 
several different fluids (author?) [53|. The two main 
interests for such a choice are the following: 

i) The effective thermodynamic exponents have been 
precisely fixed at (non Ising) values of 7 eos = 1-190, 
Asos = 0.355, a eos = 0.100 (the subscript eos recalls the 
origin of these effective values). As shown in Fig. |4j 
the value 7 eos = 1.190 is precisely within the selected 7 e 
range beyond the Ising-like preasymptotic domain, but 
well inside the extended asymptotic domain 7 e < 7 £, 
which corresponds to At* < C^\ B ; 



ii) The effective values of the thermodynamic ampli- 
tude r+ s [see below Eq. (|B8[) ] of the isothermal com- 
pressibility were then obtained only using two adjustable 
(fluid dependent) parameters (namely k and a), which 
are the two characteristic parameters involved in the 
scaled equation of state. As shown in Fig. 01 "equivalent" 
values of T+ (7 eos ) at 7 eos = 1-190 can be simultaneously 
obtained by a scale transformation between the point P 
(on the physical curve) and the point M (on the master 
curve) which also involves only two characteristic param- 
eters (namely Y c and Z c ) [admitting that the parameter 
A* e which accounts for quantum effects is known]. 

Therefore, both in quantity (two), and in nature (Ising 
like), the fluid dependent parameters k and a appear 
"equivalent" to Y c and Z c , except the noticeable distinc- 
tion of their respective determination, outside the Ising 
like preasymptotic domain for the {k; a} pair, asymptot- 
ically close to the critical point for the {Y c ; Z c } pair. 

Now we compare the respective values of T+ s and 
r"e (7eos) for twelve selected fluids. From the linear model 
of the parametric equation of state, Eq. (|B7|1 can be 
rewritten as 

^ = r+ 8 (AT*r 7M " (B8) 

where r+ s is related to the characteristic parameters k 
and a as follows 

r+ s = - a (B9) 

Considering then the restricted form of the linear model 
such as analyzed in Ref. (author?) [59|, k can be esti- 
mated from the relation 

where 6g LH = 1.3908 is an universal quantity while xq is 
a fluid-dependent parameter related to the value of the 
effective amplitude of the coexistence curve (associated 
to the value /3 eos = 0.355 of the effective exponent). The 
values of x and a can be found in Ref. (author?) [H§|. 
They are reported with the corresponding k values in 
Table [V] (columns 2 to 4, respectively) for the selected 
twelve fluids (column 1). The related values of r+ s ob- 
tained by using Eq. (|B9 j) are given in Table [V] (column 
5). 

The unequivocal scale transformation between the 
points P and M is given by the relation (see Fig. [4]) 

A* <Y\~ le °" 

r+( 7eos ) = z+ e ( 7eos ) (Bii) 

where Z+ e (7 e0 s) is the effective master amplitude for the 
T* ( 7e ,i/)-value satisfying the condition 7e ,i/ = 7eos = 
1.19. For practical use of Eq. (jB 1 1|) . the crucial advan- 
tage is given by the uniquivocal scale transformation be- 
tween the points T (on the theoretical curve) and M (on 
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Fluid 


xo 


a 


k 


r + 

eos 


Y c 


z c 


It (7eos) 


At* (7eos) 


r%{Vt) 




(author?) [59] 


(author?) [591 


Eq. (|B10P 


Eq. lfB9jl 






Eq. flBllJ 


Eq. (|B13[) 




3 He M 


0.489 


4.63 


0.9235 


0.1995 


2.3984 


0.30129 


0.20003 w 


2.326 10~ 2 


-0.283 


Ar 


0.183 


16.5 


1.309 


0.07934 


4.3288 


0.2896 


0.09284 


1.289 10" 2 


-14.5 


Kr 


0.183 


16.5 


1.309 


0.07934 


4.9437 


0.2913 


0.07887 


1.128 10~ 2 


0.6 


Xe 


0.183 


16.5 


1.309 


0.07934 


4.9137 


0.2860 


0.08084 


1.135 10" 2 


-1.86 


2 


0.183 


15.6 


1.309 


0.08392 


4.9864 


0.28797 


0.07890 


1.119 10 -2 


6.36 


N 2 


0.164 


18.2 


1.361 


0.07478 


5.3701 


0.28887 


0.07201 


1.039 10" 2 


3.84 


CH 4 


0.164 


17.0 


1.361 


0.08006 


4.9838 


0.28678 


0.07928 


1.119 10' 2 


0.99 


C2H4 


0.166 


17.5 


1.355 


0.07744 


5.3487 


0.2813 


0.07431 


1.043 10" 2 


4.22 


C0 2 


0.141 


21.8 


1.436 


0.06587 


6.0104 


0.27438 


0.06631 


0.928 10~ 2 


0.653 


NH 3 


0.109 


21.4 


1.361 


0.07353 


6.3019 


0.24294 


0.07079 


0.885 10" 2 


3.88 


H 2 


0.100 


22.3 


1.622 


0.07275 


6.8552 


0.22912 


0.0679 


0.814 10~ 2 


7.14 


D2O 


0.100 


22.3 


1.622 


0.07275 


7.0728 


0.22783 


0.0658 


0.799 10" 2 


10.6 



Table V: Two-parameter universality of the effective amplitude of the isothermal compressibility estimated from the linear 
model of a parametric equation of state and the master modification of the theoretical function. 



the master curve) illustrated in Fig. [4j That provides im- 
mediately Z+ e ( 7eos ) - Z+ e ( 7 e OS ) (eW>) 7e 
and T* ( 7eos ) = t ( 7eos ) (eW>)~\ Using Eqs. (HQ) and 



the theoretical function Xth (t) leads to the val- 
ues t ( 7eos ) = 2.392 x 1CT 4 and Z+ e ( 7eos ) = 0.456414 
related to coordinates of the points labeled T' and T" 
in Figs. UK and [5b, respectively. Using the numeri- 



3. "Universal" approximation of the logarithmic 
singularity of effective amplitudes 

Another practical application of the point to point 
transformations given in Fig. [4] can be obtained focusing 
our attention on the logarithmic singularity of any first 

derivative ( g^ p,e ) close to the Ising-like critical 



cal values of Z x and 0W> given in Table El we ob- point) for any effective amplitude power law Z+ e (t) - 



tain T* ( 7eos ) = 5.579 x 1(T 2 and Z+ e (7 eos ) = 0.15374 
Subsidiarily, in Fig. we note that the master curve 
Xqf [T*) of Eq. I|84p has a tangent curve of slope — 7eos 
at the point M of T* ( 7eos )-coordinate which corresponds 
to the effective power law 

,-1.19 



(T*) = 0.15374 (T*)~ 



(B12) 



The values of Y c and Z c for the selected fluids are re- 
ported in Table [V] (columns 6 and 7, respectively) . The 
estimated values of r+( 7eos ) using Eq. IjBlip are given 
in Table FVl (column 8). Each physical curve k* t (At*) of 
Eq. (|66|) have a tangent curve of slope — 7e0 s at the point 
of AT*-coordinate: 

At* ( 7eos ) = T^A (B13) 

1 c 

(see column 9 of Table [V| , which corresponds to the ef- 
fective power law T+ (At*) = r+( 7eos ) (At*) 7eos . 

The residuals r%(T+) = 100 ( r +^ - l) (see col- 
umn 10, Table [V]), generally lower than the typical ex- 
perimental uncertainty estimated to 10%, confirm that 
the universal features observed beyond the Ising-like 
preasymptotic domain but within the Ising-like extended 

,+ /■ + .{!/} f 

BAD — EAD — rJ 

iJ — Y a ~ EAD 

0.07-0.1, are 



asymptotic domain, i.e., At* < 



with t+ AD 



9 {l/}£+.{l/} „ nd £+,{!/} 

EAD anu i- EAD 



well-characterized by the two critical scale factors Y c and 
Z c of each fluid /. 



^-}p\ estimated from any crossover function Fp (t) given 

in Ref. (author?) Jlliwith e P , e (t) = - dLn ^ t)] ) (see 
Refs. (author?) |13l.l48l| for detail). For the suscepti- 
bility case, the logarithmic singularity of ( ff x,e ) 

extrapolated beyond the Ising-like preasymptotic domain 
is illustrated by the curves labeled ar, «m, and ap in 
Fig. IH For better evaluation within the preasymptotic 
domain, the related amplitude singularity in terms of the 
thermal field dependence is given in Fig. [5b for example 
by the curve T e of equation 



x [1 - log (t A )] t A 



(B14) 



We can approximate Eq. (|B14[) by the following "univer- 
sal" power law 



Z1 



z (z+) 



(B15) 



where Zq = 3.7 ± 0.1 and z = 0.45 ± 0.035 are inde- 
pendent of the property and the domain. The exponent 
condition z < A, leading to the ^ = 1 — u < 1, is 
conform to the logarithmic singularity of the first deriva- 
tive ( dy X 'a ) ^ ere approximated by a power law 

(X ( 7 — 7e th ) ". For practical use, we 



V^.th J 



7e,th^7 
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arbitrarily choose u = a, leading to define z = A (1 — a). 
The validity of this approximation is illustrated by the 
curve T z in Figure 03 ■ 

Correspondingly, in Figure[5ji, the physical asymptotic 
representation of T+ — T + is now approximated by the 
curve P z of asymptotic equation 



r+-r H 



Z r+a+ (At*) 2 



(7 



x! 



tf A (At*)* 



(B16) 



Using Eqs. (|B15|) and (|B16|) at f = At* = 1, we obtain 



r+ (i) - r+ = x* x $ A z+ (i) - 



(B17) 



The point to point "universal" transformation which ap- 
proximates the logarithmic singularity is then illustrated 
by the two correlated points (symbol x) at t = At* = 1, 
in Figures EJd and EJi, respectively. As expected from 
Figure HI this transformation is given by the product 

Obviously, to close the asymptotic behavior within the 
Ising-like preasymtotic domain we can also consider the 
respective asymptotic curves labeled Ta and Pa in Figs. 
UK and |St of equations 



7 - 7e,th 



= AZ^t A 



7 - 7e,exp 



= Aa+ (At) a 



(B18) 



(B19) 



Here, the point to point transformation at t — At* = 1 
(symbol x), is given by the scale factor universal power 
law i9 A . 

The above approximation of the logarithmic singular- 
ity has a practical importance for better analysis of ex- 
perimental data when the value 7 e is found in the range 



7e = 1.21 — 1.24, i.e. a value which "approaches" the 
theoretical Ising value. As a typical example we use the 
value je, P VT = 1.211 ± 0.025 obtained by Levelt-Sengers 
et al (author?) [H^| from their analysis of the pVT 
measurements of Habgood and Schneider (author?) [53| 
in the temperature range 0.2 K < T — T c < 1.8 K , 
i. e. AT*, in = 6.9 x 10- 4 , At* 1&x - 6.2 x lO" 3 and 
(At* vt ) = 2.07 x lO- 3 [see line #1, column 4, Table IV]. 
This result is then centered near to the Ising-like border- 
line of the gray paint domain previously analyzed. As 
evidenced by the matching of the corresponding points 
labeled 1 with the curves Pa and P z in Figs. \5jp and[5]l, 
such a result also appears correctly accounted for using 
the above approximation. Therefore, using Eqs. (|B19j) 
and 1|B16|) , we can easily calculate the two values of the 
true confluent and leading amplitudes of the two-term 
Wegner expansion from the following equations 



X \ P VT 



\pVT 



= 1.26666 

r+ 

1 ejVT 



1+1(7" 



At* 



(B20) 



(B21) 



= 0.057355 

These above values are in excellent agreement with the 



estimated ones 



,+ 



1.23397 and T+ 



0.057821 



(author?) [48j from application of the scale dilatation 
method. In Eq. (|B21|) . we note the practical importance 
of the prefactor Z . 

In conclusion, using Figs [4] and [5j we have explicitely 
demonstrated that the two dimensionless scale factors 
Y c and Z c , (or alternatively but equivalently dc (= i?) 
and Xq £ ), which characterize each one-component fluid 
/ belonging to {1/ }-subclass, can be used to calculate 
the isothermal compressibility over a Ising-like extended 
asymptotic domain At* < C^ AB . 
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